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Abstract 

The contribution of quarks with masses m ^> Aqcd is the only part of the structure functions in 
deep-inelastic scattering (DIS) which is not yet known at the next-to-next-to-leading order (NNLO) 
of perturbative QCD. We present improved partial NNLO results for the most important structure 
function Fi(x, Q 2 ) near the partonic threshold, in the high-energy (small-*) limit and at high scales 
Q 2 3> m 2 \ and employ these results to construct approximations for the gluon and quark coefficient 
functions which cover the full kinematic plane. The approximation uncertainties are carefully 
investigated, and found to be large only at very small values, x < 10 -3 , of the Bjorken variable. 



1 Introduction 



The production of heavy quarks in deep-inelastic lepton-hadron scattering (DIS) is an important 
process as it sheds light on a number of interesting issues in the theory of the strong interaction, 
Quantum Chromodynamics (QCD). First of all, the production mechanism can be described by 
standard renormalization-group improved perturbation theory for sufficiently large values of the 
momentum transfer Q 2 between the lepton and the hadronic states. For this purpose it is a great 
advantage that one can rely on the well-developed framework of the operator product expansion 
(OPE) which allows to understand essential features of heavy-quark DIS. In addition, the presence 
of the heavy-quark mass m introduces a second hard scale into the problem, m 2 ^> Aq CD , where 
Aqcd denotes the QCD scale. As a consequence, at any fixed order in perturbation theory (and 
neglecting bound-state effects), the heavy-quark mass defines the production threshold for the 
squared center-of-mass (CM) energy 5, e.g., S > Am 2 for heavy-quark pair production. Moreover 
the heavy-quark mass acts as a regulator of collinear divergences, giving rise to large logarithms of 
the ratio Q 2 /m 2 at Q 2 ^> m 2 , i.e., in the kinematic regime to be considered when matching QCD 
with rif light quarks and one heavy flavour to a theory with rir + 1 light quarks. 

Experimentally, heavy-quark production in DIS has been studied in fixed-target experiments 
and, in particular, at the electron-proton collider HERA. In neutral current reactions, measured with 
high accuracy at HERA, a considerable part of the inclusive DIS cross section at small Bjorken-x 
is due to the production of charm quarks, see, e.g., Refs. [HJ-|3]|. In this kinematic regime heavy- 
quark DIS is dominated by the photon-gluon fusion process y*g — > ccX, and the respective high- 
precision measurements can provide invaluable information on non-perturbative parameters in the 
cross sections such as the gluon distribution of the proton and the strong coupling constant a s . 
Given sufficiently accurate data, they can even facilitate determinations of the heavy-quark masses 
in a reaction with space-like momentum transfer. With more data analyses from HERA run II 
still to be finalized, the constraints from heavy-quark DIS on global fits of parton distribution 
functions (PDFs) and its resulting consequences for collider phenomenology are, perhaps, the 
most important aspects in the era of the LHC. 

It is thus of great importance to provide as accurate theoretical predictions for heavy-quark DIS 
as possible. Within the standard perturbative approach, this requires consideration of higher-order 
radiative corrections. At present the theory predictions for neutral current heavy-quark production 
rely on the complete next-to-leading order (NLO) QCD corrections [4] which are often used via 
the parametrizations of Ref. [0, see also Ref. [|6) for a check and minor corrections. The complete 
contribution of the next-to-next-to-leading order (NNLO) is not known. However, partial informa- 
tion in various kinematic limits has been accumulated over the years. This includes in particular 
the logarithmically enhanced terms near threshold which are accessible by means of the soft-gluon 
exponentiation Q. Also the high-energy limit of heavy-quark DIS has been known from small-x 
resummation to all orders at leading-logarithmic accuracy for a long time 0. In the asymptotic 
regime of Q 2 3> m 2 fully analytic results have been obtained, and the NLO calculations ll9lfT0l have 
been extended recently towards NNLO by the computation of a number of lowest even-integer 
Mellin moments [fTTT - TOl . see also Ref. [|T4||. Finally, the dependence on the renormalization and 
mass-factorization scales is known exactly at NNLO from standard renormalization-group argu- 
ments Elm and can be used as an independent consistency check in all these limits. 



1 



In the present article, we provide approximate NNLO QCD corrections for the heavy-quark 
part of the dominant structure function F 2 in photon-exchange DIS. This is achieved by extending 
and, for the first time, combining the available NNLO information from the different kinematic 
regimes, i.e., from threshold (s ~ 4m 2 , where s is the partonic CM energy), high energy (s ^> Am 2 ) 
and large scales Q 2 ^> m 2 , in particular for the dominant gluon-initiated contribution. Specifically, 
we present results for all four soft-gluon enhanced logarithms near threshold (see Ref. |[T6l for 
a previous brief account), and we employ the all order-result in the high-5 limit [8] to derive an 
analytic expression for the leading In s term at NNLO. Finally, the known Mellin moments of the 
heavy-quark operator matrix elements (OMEs) |[TTI - fT3l are employed together with the three-loop 
results for massless quarks [fTTI to construct approximate x-space expressions for the heavy-quark 
coefficient functions at Q 2 ^> m 2 . By combining these individual constraints we construct NNLO 
coefficient functions for heavy-quark DIS which, while still being approximate, represent the most 
comprehensive results possible at this point. Below we will provide a detailed documentation of 
the required calculations as well as estimates of the accuracy of these approximate results. The 
resulting improvement of the predictions for heavy-quark DIS and the low-<2 2 small-* limitations 
of the present NNLO results are then illustrated in a brief phenomenological study. 

The remainder of this article is organized as follows: In Section [2] we introduce our basic no- 
tations. Section [3] is devoted to the calculation of the new NNLO results for the gluon coefficient 
function near threshold (Section I3TTT ). at small-* (Section 13^21) and for scales Q 2 ^>m 2 (Section l33l) . 
In the latter two cases we also address the light-quark coefficient function. In Section |4] we com- 
bine these results and provide an approximate NNLO expression for both coefficient functions for 
which we perform a number of quality checks. Section |5] addresses their impact on the charm- 
and bottom-production structure functions F 2 C and F 2 b at NNLO and the limitations of the present 
results. Finally we conclude and briefly discuss possible future improvements in Section |6] The 
appendices contain, in Appendix lAl the derivation of the analytical result for the gluon coefficient 
function at NLO in the limit s ~ Am 2 and, in Appendix|Bl the (up to two OMEs) exact, if somewhat 
lengthy results for the heavy-quark coefficient functions at asymptotic values Q 2 ^> m 2 . 



2 Setting the stage 

We start by briefly reviewing the main results on heavy-quark production in deep-inelastic scatter- 
ing mediated by neutral-current exchange. For the scattering of a charged lepton e off a proton P 
this reaction, 

e(l)+P(p) e(l')+q h ( Pl )+q h (p 2 )+X, (2.1) 

is dominated by virtual photon exchange if the momentum transfer Q 2 = —q 2 = — (/ — I') 2 is 
much smaller than the Z-boson mass, Q 2 <C M|. The quark pair q h q h in the final state is heavy, 
m 2 ^> Aq CD . The cross section, integrated over the momenta of the outgoing heavy (anti-) quarks, 
is written in terms of the heavy-flavour structure functions Fk(x, Q 2 , m 2 ) with k = 2,L as 




(2.2) 
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where a is the electromagnetic coupling constant, m the mass of the heavy quark, and the well- 
known DIS variables x and y are defined as x = Q 2 / (2p ■ q) and y = (p ■ q) /(p ■ I), respectively. 

Disregarding power corrections, the heavy-quark contribution (12.21) to the DIS structure func- 
tions can be written in terms of a convolution of PDFs and coefficient functions as, see Ref. 0, 

2 z max A f \ 

F * (x ' GW) = W ^ / T/ifp^j^hW^^), (2-3) 

47T 171 i = q,q,g J * Z \Z J 

where z max = 1/(1 +4m 2 /<2 2 ) and e h is the heavy-quark charge. The strong coupling constant at 
the renormalization scale /j r is denoted by a s , and the PDFs for the parton of flavour i are f^x./u 2 ) 
at the factorization scale /Uj. The kinematic variables T| and £, in Eq. ( 12.31) are given by 

n = A" 1 > S = % > (2-4) 

and the partonic CM energy is s = Q 2 (l/z— 1). Instead of T|, which measures the distance to the 
partonic threshold, one often uses the heavy-quark velocity (3 or the variable p, 



(3 = ^l-Am 2 /s , p = Am 2 /s . (2.5) 
The coefficient functions of the hard partonic scattering process can be expanded in a s as 

co oo J 2 

CkM&f) = I ^o^fo,^ 2 ) = £ (47ta.y J>#*°(ti,$) In'*-, > (2-6) 

7=0 7=0 /=0 m 

where we have identified the renormalization and factorization scales, /j = /j^ = /j r . This can 
easily be undone by expanding 0C s (/jy) in terms of CL S (/J r ) using the standard QCD beta-function. 

In the above normalization, the coefficient functions at the leading order (LO) read 



4%,Q = 4%X)+c^( % i), (2.7) 

c S ( " ,e = 2(1+^+4/4)3 {-2(d+n-?/4)^+n + i)P 

+ (2(l+l 1 ) 2 +5 2 /8 + 2ri + l)L(P)}, (2.8) 

«Sd.e = 2(1+ T+W { 2p( ' +11> " L(P) } <2 ' 9) 

with L((3) = ln((l + P)/(l — P)). Here Cj g and c|°j denote the contributions from transverse 

and longitudinal photon polarizations, and 7> = A for the colour group SU(N C ). The radiative 
corrections to Eq. (12.61 ) at NLO have been known for a long time from Ref. fl4]. Unlike the massless 

coefficient functions, the massive NLO coefficient functions c^j cannot be expressed in a simple 
analytic form. Instead, compact parametrizations for these functions were presented in Ref. [0, 
with minor corrections provided later in Ref. [|6] (see Ref. [|20l for parametrizations in Mellin-A^ 
space including complex values of N as required for a numerical Mellin inversion). 
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The NNLO coefficient functions cl •' in Eq. (12.61 ) are not fully known, although, as outlined 
above, substantial partial information is available. On the other hand, all scale-dependent terms at 

(2 1) (2 2) 

this order, i.e., c k - and c k - , have been constructed by means of renormalization-group argu- 
ments B71I2TJ and are completely known in numerical form, using the parametrized NLO results 
and the well-known expressions for the NLO splitting functions, see Refs. ir71[T5l. In the following 
we will thus focus on improved predictions for the scale-independent parts of the NNLO coeffi- 
cient functions. To be precise, we will confine ourselves to the gluon coefficient function 
which is by far the most important contribution to the heavy-quark structure function F2 after the 
convolution with the gluon distribution in Eq. (12.31 ). and to the pure-singlet light-quark coefficient 
function cl 2 ' " 1 for the heavy-quark contribution proportional to e 2 . 

For further reference we finally introduce the short-hand notations 

Ln = ln , Lr, = ln ^ and a s = , (2.10) 
m l fj 4% 

and we also note that throughout this article we will employ the on-shell scheme for the heavy- 
quark mass m, i.e., the pole mass. Heavy-quark DIS with a running mass in the MS scheme has 
been considered in Ref. lfT5ll . 



3 The NNLO corrections in three kinematic regions 

In this section we provide improved predictions for the above two NNLO coefficient functions 
c^g ^ and c"2q \ We address, in this order, the threshold region s ~ 4m 2 , the high-energy regime 
s S> 4m 2 , and the high-scale region Q 2 ^> nr. In all these limits the higher-order corrections 
exhibit well-known features and regularities which we briefly discuss and exploit. The results will 
be assembled into approximate coefficient functions for the whole kinematic plane in Section 31 



3.1 Threshold limit and soft-gluon resummation 

Partonic cross sections generally receive large logarithmic corrections near threshold, which appear 
as the result of an incomplete cancellation of the soft (and collinear) contributions between the 
real and virtual corrections due to the reduced phase-space near threshold. In the present case of 
massive quarks, the corrections appear as a function of the heavy-quark velocity (3 at each order of 
the perturbation series in the form a" ln ra [3 with m < In. The production threshold corresponds to 
13—7-0, recall Eq. (12.51 ), and the highest logarithms can be resummed to all orders in perturbation 
theory, see, e.g., Refs. [|22] - l26l . In addition to the logarithmic terms, heavy-quark pair production 
also receives corrections of the form a"/(3 _m h/|3 with m < n at all orders from the Coulomb- 
exchange of gluons between the heavy quarks. Also these terms can be resummed fl27l . 

The raison d'etre of phenomenology based on threshold resummation derives from the well- 
known fact, see, e.g., Refs. [71|28l|29]|, that, at not t0 ° large values of Q 2 , the convolution of 
the coefficient function for F2 with the gluon PDF in Eq. (12.31) is often dominated by rather low 
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partonic CM energies. Therefore the NNLO predictions of the threshold resummation can provide 
useful information on this numerically important contribution to F 2 . Previous research has already 
determined exactly the two highest Q and approximately the third [1301 threshold logarithms at 
this and all higher orders. In what follows, we derive exact expressions for all four logarithmically 
enhanced terms together with the complete Coulomb corrections for the NNLO gluon coefficient 
function cj_ ' in Eq. (12.61) . a result that we have already reported in a numerical form in Ref. IfToll . 
The (light-) quark coefficient function c 2 q for heavy-quark DIS is non- vanishing only from NLO 
and does not exhibit any enhancement near threshold. 



According to the general formula for the threshold resummation 11241 - 12611 . the gluon coefficient 
function in Mellin space (recall Eqs. (12.41) and (12.51) for the definitions of ^ and p) is of the form 

c 2jg (a s ,A0 = jTdpp^c^&P) 

= 4° g V) • So(«s,A0 ■ exp [G N (a„m 2 )} ■ (l + O^ln^)) . (3.1) 

The coefficient cffg(N) is the Mellin transform of the LO coefficient function of Eq. (|2.71) . In the 
large-/V limit (which is dominated by the transverse component c T ) it is given by 

c S (Ar) = (TT^4) TtPr- (' + ■ (3.2) 

g (a s ,N) denotes the matching function to be discussed below, cf. Eq. (13.71) . The exponential 
factor resums the threshold logarithms to all orders in oc s . The exponent can be written in the 
standard form [TEH 



rl 7 N-1 _ i ( r 4m 2 (l-z) 2 A a ^ 

Gn = J dz ~T^z~ 1 J 2 ~lf A z ^ q2) "> +D r*-*»* ( as ( 4m t 1 " z ] )) } • (33) 

Here the first term includes the gluonic cusp anomalous dimension A g , known to order B2U33H . 
which governs the process-independent soft-collinear gluon emission off the initial gluon. The 
second term, Dy* g ^ qh q h , collects process-dependent effects of soft-gluon emission from the initial- 
and final-state particles and is built as 

D rg-+4hqh( a s) = i£g(as)+A^(a s ) (3.4) 

from D g and D qhq - h which are, respectively, the soft anomalous dimension for Higgs production in 
gluon-gluon fusion (known to order a s 3 ll34l[35l ) and for the heavy-quark production in the color- 
octet channel (known to order a s 2 [|36] - [38l . see also Refs. |[39l|40]l). For reference we assemble 
Dy*g^ qh q h , employing the usual notation (12.101) for the strong coupling, 

D Tg^q h q h ( a s) = -4C A a s + 4C A C A ^-^ + ^^ 2 + 5^ 3 ) +n f (57-5^2) a s 2 + 0(a s 3 ). 

(3.5) 

Here and below L, n stands for values of the Riemann zeta-function. 



5 



To next-to-next-to-leading logarithmic (NNLL) accuracy, as needed for the NNLO corrections 
we are considering, the exponent in Eq. (13.31) takes the form 



\nN 8l (X)+g 2 (X)+a s g 3 (k) + ... (3.6) 



with X = $oa s \nN and, again, a s = a s /(4n). The functions gj (i = 1,2,3,...) are defined such 
that gf(0) = 0, with the first term InN g l resumming the leading logarithms and so on. Explicit 
expressions for gi can be obtained from Ref. OTI . even with the /j r and /j^ dependence separated 
(for the computation see, e.g., Refs. H41U42II ). with the obvious replacements A q — > jA g , D q ^^D g 
and Dqq — > —D qh q h in Eq. (A. 5) - (A.9) of Ref. OTTl (note especially the unconventional sign 
for the soft anomalous dimension of color-octet heavy-quark production, there denoted by Dqq). 
Below we will expand in terms of ln/V = lnAf + ln y e , which absorbes the otherwise ubiquitous 
Euler-Mascheroni constant into the logarithms. 

The process-specific new information needed for heavy-quark DIS resides entirely in the match- 
ing function g (oc s ,/V) in Eq. (13.11) which we need to discuss next. This function can be expressed 
as a product of a hard coefficient g[)(oc s ) and a Coulomb coefficient g c (a Sl N), where the latter 
specifically induces a dependence on the Mellin variable ./V, 

g (a s ,AO = g£(a s ) g%(a s ,N) 

r, h(l) , 2 h(2) , \ /1 c(l) . 2 c(2) , x 

= (l + a s g { ' + <g v + ...) (I +a sgQ { ' + < g Q { ' + ...) 

= l + a s (g h U + g^) + al(g^ + g^^ (3.7) 

Such a factorized form for the Coulomb corrections and the threshold logarithms has been dis- 
cussed for top-quark pair production at hadron colliders before H43II441 . see also Refs. [|36ll431l for 
studies in the framework of soft-collinear effective theory. For heavy-quark DIS a similar factoriza- 
tion applies to the gluon coefficient function, since the structure of radiative corrections to the DIS 
subprocess y+ g — > q\ x + q\ x + X is similar to top-quark hadro-production in g + g — > q/ r + q~h + X, 
i.e., essentially the same up to replacements of the corresponding color factors. 

The matching coefficients are fixed by comparing the (Mellin- inverted) expansion of Eq. (13.11 ) 
to known fixed-order results, in our case those of Ref. [@]. The calculation required in this case is 
discussed in Appendix lAl and leads to 

g c W (N) = 2^[C F -^) X I^-, (3.8) 



g Ml) = (41n 2 2-41n2-40 + 27r 2 )C A +c ^) + [8C A -c (^)]^ , (3.9) 




where 



c ® = C A {50-^ + 12 ln ^^ ^ / ^ +41n 2 ( v ^(v-l)/2)+ln 2 (l+V2) 

+ 61n(2 + V2)-41n 2 (2 + ^/2)+2Li 2 (-^) + ^-4i^±|^ 
64 M2+j/2) _ 1 2g ln(2 + $/2) _ 160 _ ^ ln(2 + $/2) 



+ (2 + ^) 2 (4 + ^) (2 + $)(4 + $) (2 + ^(4 + ^) 
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+ 



128 



+ C P 



(2 + 5)(4 + 5) 2 
2 



-12 



4+C2 
4 + 5 



8- 



In' 



4+5 



+ 



64 



(4 + 5) 



(3.10) 



18--7T-24- 



ln(^(j-l)/2) 



81^(v / ^-l)/2)+21n 2 (l+5/2) 



-61n(2 + 5/2)+4Li 2 (- 
ln(2 + 5/2) 



+ 128 



2 + 5 ; 2 + 5 
544 



48 | g hi(2 + 5/2) | 



360 



(2 + 5)(4 + 5) (2 + 5)(4 + 5) 2 



+ 48 



ln- 



2+5 (2+5X4+5) 

5(y-l)/2) ln 2 (l+5/2) 



4+5 



+ 



44 + 40^2 

4+5 



120 ln g + ^+256- 



(5) 



4C A 



2 + ln 1 + 



(2 + 5) 2 
4 



(2 + 5) 2 (4 + 5) 



16 



8 

(4 + 5) 



4 + 5 



272 



(4 + 5) 2 

(3.11) 



with y = -y/l +4/5 and C A = 3, C F = 4/3 in QCD. The additive constants taken out of c (5) and 
c (5) in Eq. (13.91) are conventional; the reason for the above choice will become clear below. 

The Coulomb coefficient gQ(a Sl N) in Eq. (13.71) can be obtained to order a 2 by making use 
of the result of the non-relativistic (NR) cross section in e + e~ — > q h q h calculated up to this order 
in Refs. H6ll47l . The only difference in the present case is the color structure, i.e., the Coulomb 
corrections for the color-octet state require the colour factor replacement C F — > (C F — C A /2) in 
the corresponding e + e~ result [|46ll47l . Thus we identify the NR part of the NNLO coefficient 
function as 



,(2) 
%8 



(5,P) 



NR 



.(0) 
'2 ; g 



(5,P) 



(An) 4 



C 



c A 



? C A-?"/- 4 Poln(2|3)-2po^ 



32C F 7t 2 lnp 



+ o(P) 



(3.12) 



where |3o = 11/3^ — 2/3n^ is the leading coefficient of the QCD beta function and we have 
removed the term corresponding to the product of the NLO hard and NLO Coulomb corrections 
to avoid double counting. Note that the single logarithmic term originates in the non-leading part 
of the NR QCD potential ||48l and is not included in the soft-gluon exponent Gjy. Performing the 
Mellin transform, we obtain ™ as 



c(2) 
So 



TC* C 



2 

IAN 



/62 20 \ AN 

[jC A --n f -A^\n2-2^L^ ^j — 



]^-\nN + Uc F - ^ J n 2 - (2N) -32C F ( -UnN+l -ln2 



.(3.13) 



The only missing information for the complete NNLL + NNLO resummation is the coefficient 
g Q in Eq. (13.71) . which requires the presently unknown NNLO constants in p\ i.e., the two-loop 
analogues of c (5) and c (5) in Eqs. (13.101) and (13.111) , which can be determined only from a full 
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three-loop calculation of the heavy-quark contribution to DIS in the soft-gluon limit. Setting this 
constant to zero in [3-space, we find 



* (2) 



135416 
27 

932 9 



In 2 



35888 , 2 , 11672 Y 
9 In 2 H — C,3 

488 7 .2 



88856 . 1340 r . 16 4 
- H 9 — 71 + — 7t 



27 



21584 



+f 7t 2 In 2 - 1 384 ^3 In 2 + 7t 2 In 2 2 + In j 2 - 640 In 4 2 J C| 

Z4592 224 r 80 Jl 6512 416 ,3 , a ^.2 ,„ , 1376 ,2 V „ 

+ I 27 — "n - C,3 — -Q- 7C yf In 2 — -g- In 2 + 8 7t In 2 H — g— In Z\C A n 



+ (-32 In 2 2 + 56 In 2 + 2 7t 2 - 40) C A c (£) 

- "(-^7t 2 ln2+ 2881n 3 2+^ln2+ 347t 2 - ±^ ln 2 2+ 224^ 3 - ^) C 
+ (320 + 64 ln 2 2 _ 416 ^ _ 4 ^ Q ^ + (8 ln2 - 8) C A c (^) 



+ (2tc 2 - 32 ln 2 2+ 56 ln2- 40) C A c (5) L p + (An 2 - 32 ln 2 2 - ^ + ^ ln2^ C 2 
(f - f ln2) C A n f + (8 ln2 - 8)C A c (5)] L 2 (3.14) 



+ 



which cancels all NNLO constant terms originating in the inverse Mellin transform of the iV-space 
resummation formula (I3.ll ). 

Supplemented by an estimate of the missing [3-space constant, e.g., from a Pade estimate as in 
Eq. (13.191) below, the results in Eqs. (13.61 ) - (13.141) can be used to provide an approximate NNLL 
+ NNLO exponentiation of the gluon coefficient function in heavy-quark DIS, which is analogous 
to the case of hadronic top-quark pair production in Refs. H26LI3 111 . 



In the present article, though, we are rather concerned with improved predictions for c 2 g to 
NNLO, essentially using Eq. (13.11) to generate results to fixed order in perturbation theory. To that 
end we now write down the (3-space threshold behaviour for c 2 g up to NNLO. The LO and NLO 



results 



are 



JO) 

'2,g 



&P) 



P 



1+^/4 



+ o(p 3 



(3.15) 



and 

C 2,g 

with 



5,P) 



JO) 

'2,g' 



S,P) 



(4tc)2 

_ (l)const 



1 6 C A In 2 (3 + [48 C A In 2 - 40 C A ] In (3 + ( 2C F - C A ) ^- + 8 C A In p L^, | 



+ c£ M "~&P) + 0(P*) 



(3.16) 



Jl)const 

'2,g 1 



S,P) 



JO), 
'2,g 



4,P) 



(4ny 



|c (^) +36C A ln 2 2 - 60C A ln2 + L^[8C A ln2-c (^)]} (3.17) 



in terms of c (^) and c (i;) in Eqs. ( 13.101) and (13.1 II) . The convention for these constants derives 
from Refs. BUO, where the logarithmic enhancement was expressed in terms of In (8|3 2 ), £ = 1,2, 
and ln(4J3 2 ) for the scale-independent and scale-dependent parts, respectively. 
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(2) 

The corresponding NNLO result g is obtained from our results above using the Mellin trans- 
form in Appendix A of Ref. OTTl . 



J2) 
'2,8 



&P) 



JO) 

'2,g 



S,P) 



+ 
+ 



(4tc) z 



768 ln2 



128Ciln 4 (3 



(3.18) 



6464 \ 
"9" J 



128 



(l7281n 2 2-3232 ln2- 



( C A n 7 + 128^] ln J B 

208 , 15520 \ 
H 



r 2 



+ (64 ln2 - ^ ) C A n f + 16 C A c (5) + 32 C A (c F 



P 



5121n2 + 



1136 



32 



15520 



ln2-208^ 2 ln2 + 936^ 3 + — n 



(l7281n 3 2-48481n 2 2 + 

88856 \ . / , 2. 640 
-^7-) C ^ + ( 961n2 -- 

+ (481n2-40)C A c (^) + {(-y + 32rn2)c A + ^n / }(c F -^)^ 

104 



In 2 



16 2 4592 



C A n f - 



6721n 2 2 + 9761n2 + 



32C>(C 
C, 



C, 



7T 



4160 



■Jt - 



Q 2 + 



321n2 + 



320 



C A n f 



+(48 ln2-40)C x cb® - 8C A c (^) - 16C A (c F - ^) ^ }l„ 



+ <^ 641n2 



s -(5)}^ 2 ]lnP 



Q + -Qn / -8C A c ( 



~ / n 68 8 o 

+ ^81n 2 2-y ln2 + -7l 2 
/26 
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658 



C A + (1^2-^)^ + 200^) 



+ 



2co(^) L 



c. 



% 2 4 
T + 3 



c 



27^ 



+ o(P). 



Several checks can be performed of this result which provides, for the first time, all logarithmically 
and Coulomb-enhanced contributions. First of all, the scale dependence at NNLO is fully known 
lP71[l~5Tl from the renormalization group (see also Section |4j, and we have verified that all terms 
proportional to powers of in the threshold expansion in Eq. (13.181 ) agree. Another strong check 
is the agreement of Eq. (13.181 ) with the results presented in Eq. (A.l) of Ref. Il38l for the hadronic 
heavy-quark production in the gluon-gluon-fusion channel, gg — > q^qh, after replacing one gluon 
by a photon (together with the appropriate substitutions of the colour factors). 

The NNLO constants in B multiplying the Born term are currently unknown. This includes the 
scale-dependent term Cj g T|) which is only known in a numerical form beyond the threshold 
logarithms. Based on Eq. (13.171 ), one can write down a [0/1] Pade estimate, 



c (0) (£ B) 

- H^fT ( c o(^) + 36Q In 2 2 - 60C A ln2 + L„ [8C A ln2 - c © 



_ (2)const 



S,P) 



(3.19) 



of which we will employ the scale-independent part at low £ in Section [4] below. 
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10 4 10 3 10 2 10 1 1 10 10 4 10 3 10 2 10 1 1 10 io 2 

r| r| 

Figure 1 : The successive logarithmic threshold approximations for the NNLO gluon coefficient function 
in Eq. ( 13.181 ) and the 1/(3 Coulomb contributions at two typical values of £, = Q 2 /ni 2 depending on r| 
defined in Eq. (12.4b - The two highest logarithms were already known from Ref. Q. 

The main new NNLO result of this section, the scale-independent part of Eq. (13.181) . is illus- 
trated in Fig. Q] for two values of Q 2 . It is clear that a stable low-T) result is established only by 
including our new results for the third and fourth logarithms. As a first estimate of the remaining 
uncertainly at r\ > 1 for low values ^ > 1 we have included the Pade estimate (13.191) , assigning a 
100% uncertainty (its small value at ^ = 20 is accidental and irrelevant in what follows). 

3.2 High-energy expansion 

Let us next discuss the high-energy limit of the coefficient functions c 2 g and c 2q . To that end 
we build on the results of Ref. (HI from the small-* resummation, which had been derived even 
before NLO corrections were available and thus served as a check on Ref. 01. At fixed Q 2 the 
high-energy limit s = Q 2 (\/x- 1) -»■ °° implies x -»■ 0, thus r\ ^/(4x) + 0(1) due to Eq. (|24|) . 
Therefore we consider the coefficient function c 2 g as a function of x throughout this section. In 
Mellin iV-space, cf. Eq. (13.11) for the definition, the dominant contributions in this region behave as 
a"/N n and can be resummed to all orders. 

One can define the partonic cross section 6 2N in Mellin space corresponding to the gluon 
coefficient functions c 2 g (see Eq. (13.11) ) as 

m 2 6 2N = aa s e 2 q c 2:g (a s ,N) , (3.20) 
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which contributes to the structure function F 2 via Eqs. (12.21 ) and (12.3b ■ In complete analogy, if 
multiplied with the Mellin moments of the gluon PDF f g pj, g 2 n also defines the gluonic con- 
tribution to the total hadronic cross section. It is an interesting observation of Ref. [8J that the 
<2 2 -dependence in this quantity can be factorized if one considers the photo-production limit, i.e., 
the process yg — > q h q h X with a real photon. This leads to the ansatz 

°2,;v(£W) = K 2 ^)o ygjN (m 2 ) . (3.21) 

The function K2,n summarizes all ^-dependence originating from the momentum transfer Q 2 of the 
off-shell photon in DIS electro-production of heavy quarks, and Oy gj N denotes the Mellin moments 
of the gluonic contribution to the total hadronic cross section for photo-production of heavy quarks, 
see, e.g., Ref. EHEOl. 

The resummation of the logarithms at high energy is based on the framework of un-integrated 
PDFs in transverse momentum k t and the concept of ^-factorization, a procedure which involves 
two steps: first computing amplitudes with the initial particles off-shell in k t , and second perform- 
ing the convolution with a gluon PDF which has the small-* corrections included. For heavy-quark 
photo-production, this allows to express Oy g ,N in Eq. (13.211) as the product (8) (cf. also Ref. 115010 . 

m 2 a Ygi jv(m 2 ) = fg,N(m 2 )h N (y N ) , (3.22) 

that is, in terms of the gluon PDF and the impact factor h N (yiy) depending on the anomalous 
dimension y N . We are interested in the perturbative regime with y N given by 

a\ a / a \ 4 



with a = C A a s /n. This well-known result for y N arises as the solution to the following condition 

on%(y), 

in a perturbative expansion for y^C 1 and requiring y(0) = 0. Here %(y) is the Mellin transform 
(with respect to the transverse momentum kf) of the lowest order BFKL kernel expressed through 
standard \|/-functions as 

X(y) = 2x|/(l)- ¥ (y)-i|/(l-y) = I + 2^ 3 y 2 + 0(y 4 ) . (3.25) 

The anomalous dimension y N governs the high-energy behavior of the gluon PDF in Eq. (13.221) as 

f g , N (m 2 ) = l-p) f g . N (n 2 ) , (3.26) 
which up to terms relevant at NNLO simply gives 

(y) =l+y N L,+^y 2 L 2 + 0(y 3 N ). (3.27) 
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In the perturbative regime the impact factor is needed only in the limit N <^y N , and can thus 
be approximated as 

MY*) = K1n)(1 + 0(N)) (3.28) 

because h N (jN) is free of singularities and terms of order N are subleading. The calculation of the 
impact factor h is performed by considering off-shell amplitudes in a fc r -factorization scheme. For 
photo-production of heavy quarks in yg — > q h q h , and with the incoming photon and gluon being 
off-shell by an amount k t , this results in |j8~l l50~ll5TI 

h(y N ) = l a a s ef^^-B(l-y N A-y N )B(l+y Nl l-y N ) (3.29) 

in terms of standard Euler Beta-functions. For y N <C 1 this expression leads to 

h(y N ) = Ttaa s e 2 q {l + ^y N + ^y* + 0(yi)} . (3.30) 

Note that the precise definition of the fc r factorization scheme for the impact factor and the conver- 
sion to the MS scheme affects the result in the perturbative regime only at order y^, see Ref. Il50l . 
This is beyond the NNLO accuracy in a s we are aiming at. 

Finally the ratio K 2 ,n between the cross sections for DIS heavy quark electro-production and 
photo-production in Eq. (13.211) has been computed in closed analytical form in Ref. flU, 

^® ^ H) V5j(l + 2 V ) {f( 1+ ^ 

+ (l+|) Y " (2 + 3Kv-3Kj-|(l+Tav))2i :, i(l-Yi„j;l;^4)}, (3.31) 

where 2 F X denotes the hypergeometric function which can be conveniently expanded for small y N 
in terms of harmonic polylogarithms H^(jc) 11521 (see e.g., Refs. II53H561I and Appendix|B]). To that 
end, we have used the HypExp2 package (571 in Mathematica, which gives us 

2z 2 F l (\-y N ^;l;z 2 ) = L(z) +y N {H(-+;z) +L(z) ln(l -z 2 )} - \y 2 N {#(-+, -;z) 

-H(-,+;z) ln(l -z 2 ) -L(z) ln 2 (l -z 2 )) + 0(y 3 ) , (3.32) 

where L(z) = ln((l +z)/(l -z)) andH(-,+;z),H(-,+,-;z) are short-hand notations (cf. 11551 ) 
for combinations of harmonic polylogarithms 

H(-,+;z) = ffi,i(z)+Hi_i(z)-H_i,i(z)-H-i,-i(z), (3.33) 
H(-,+, -;z) = Hi,i i i(z)-Hi,i_i(z)+Hi_i ) i(z)-Hi_i_i(z) 

- #-1,1,1 (z) + -i(z) " H-i _i,i (z) + _i _i(z) . (3.34) 
For the expansion coefficients of Eq. (13.321 ) we define the functions 

m) = | V / ^+4)//(-,+; V /^+4)), (3.35) 
m = y^/^ + 4)L^/^ + 4)), (3.36) 
= | V / ^/(^ + 4)//(-,+,-; V / y(^ + 4)), (3.37) 
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of which /(£) and 7(^) have been introduced already in Ref. [8] where also an expression for /(^) 
in terms of normal (di-)logarithms can be found. 

At this stage, all that remains is to perform the expansion of Eqs. (13.221) and (13.311) to second 
order for small by assembling all formulae provided so far. Converting from iV-space back to 
jc-space we obtain for the gluon coefficient function, 

4!]m = c Aii{^+(i-|)^)+(f -i) / ^)+ L 4t + ( 1 -?) / ^]} + °W' 

(3.38) 

<®M) = Cilnx3^{-i|-i(l-|)/(^)rn(l + |)-i(l3-f)/(^) 



^(71-|)7(^ + I(l-i)^) + L,[-|-f(l-i)/(^ 



(l-l) + 0(*°) 



(3.39) 



where Eq. (13.381 ) agrees with the expression given in Ref. []5] (cf. Eqs. (19) - (22) there), while the 

(2) 

result for the NNLO quantity c\ is a new result. 



For further reference (in Section 13 .3 J where we will also provide graphs of the above results), 
we have finally derived explicit results for the \ 3> 1 limits, i.e., Q 2 ^> in 2 , of Eqs. (13.381 ) and 
(I3.39I ). We obtain, at leading-logarithmic small-x accuracy, 



c A 



167T 



lL 2 + ^L Q+ f-f^+(fL Q+ l)L, +0(0, (3-40) 



104 



40 16; 



16 



" C A 2567? 



32 r 3 , 208 r2 , [2212 64 r > , 
T L Q + ~9~ L Q + \ lJ 3"'= 2 / L Q 



+ 



1472 416 
27 



128 



/ 32 T 2 . 416 j 
-5-4,2 + ^-^3+ I -T L Q + ~Q~ L Q +~Q~ ~ T^lh 



160 64. 



3 



32 



16 



+ 0(0 



(3.41) 



Here the constant terms at NNLO, i.e., the terms independent of Lg and L p (recall Eq. (12.101 ) 
above) in Eq. (13.411) are new. All other terms in this result agree with independent calculations. 
The scale-dependent L M terms can, of course, be obtained by the standard renormalization- group 
behaviour, cf. Section @] More importantly, we are also able to compare all Lq terms with the 
results obtained in the asymptotic limit Q 2 ^> m 2 below, cf. Section [3731 All powers of Lq and 
match exactly with those results, which constitutes a strong check. Another test of Eq. (13.391 ) is 
possible by comparing with the small-x limit for c 2 g derived in Ref. [|59l for the kinematic region 
Q 2 ~ m 2 (see also Eqs. (10) and (11) in Ref. Il30l ). and we find agreement at the level of a few 
parts in a thousand for Q 2 /in 2 = 1. 

In concluding this section, we briefly comment on the quark coefficient function for heavy- 
quark DIS when the photon couples to the heavy quark, that is, c 2 q in Eq. (12.31) proportional to e\ 
and the quark-singlet PDFs. Its leading small-.* term is related to that of c 2 g in Eqs. (13.381) - (|3.41l) 
by a simple colour factor replacement, c 2 q = C F jC k c 2 g , as pointed out already in Refs. Il8ll59l. 
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3.3 High-scale limit Q 2 ^>m 2 



The results in the limit Q 2 ^> m 2 build on an exact factorization of the heavy-quark coefficient func- 
tions into the respective coefficient functions with massless quarks and heavy-quark operator ma- 
trix elements (OMEs). The variable r| in Eq. (1241) factorizes as r| ->■ Q 2 /(4m 2 )(l/x- 1) + 0(1), 
and, to denote this limit, we will be considering the coefficient functions as a function of x and 
the ratio Q 2 /m 2 throughout this section. At LO, one can quickly verify for c 2g with the help of 
Eq. (12.71) that this limit reads 

2 

c^ g {x,Q 2 /m 2 ) = kx^{2(\-2x + 2x 2 )L q 



+ 6-2(1 -2x + 2x 2 )(4 + U (x)+lli(x))} + o(^j , (3.42) 



where the term multiplying Lq is proportional to the lowest-order splitting function Pqg , while 
the Lq independent term (here written in terms of harmonic poly logarithms H,n(x) ||52l , see Ap- 
pendix[B]) is proportional to the one-loop gluon coefficient function with massless quarks. 

Formally the underlying factorization may be expressed through massive partonic OMEs Aij 
(i.e., with massive quarks) and the DIS coefficient functions with massless quarks C^ lght together 
with their exact scale dependence from renormalization group invariance, see Ref. 112X11 . i.e., at 
scales [x^z Q and ^ m, 



m 2 \ light / Q 



c k MM cv (x,^,-^)= Ajt C^f ( ^ ) (x) + O ( ^ ) ■ (3.43) 



/j 2 J k >j V fJ 2 



2 

m 



Q 



Here tg) denotes the standard Mellin convolution, cf. Eq. (12.31) . The dependence on the scale /j has 
been made explicit together with the separation of scales between the OMEs and the coeffi- 
cient functions C^ ht , which is realized through \nQ 2 /ij 2 = Lq + L^. At NLO the formalism in 
Eq. (13.431 ) has first been applied in Ref. [0 to compute the corresponding two-loop OMEs for (un- 
polarized) heavy-quark DIS at asymptotic values Q 2 S> m 2 , see also the re-calculation in Ref. IfTOl . 
In the wider context, Eq. (13.431 ) implements the matching conditions for QCD with one massive 
and Hjr light quarks to QCD with «y + 1 light quarks (note that denotes the number of massless 
quarks throughout this article). The matching of course also affects the strong coupling oc s , for 
which the decoupling formulae are well known [|60ll6D . Unless indicated otherwise, we are al- 
ways working in the scheme with a s (n^), i.e., = 3 for charm quarks. Moreover, on the basis of 
Eq. (13.431) . also the all-order resummation of the logarithms in Q 2 /m 2 has been discussed in the 
literature and has led to definitions of a so-called variable flavour number scheme (VFNS) [|62~1 . 
see Refs. Il63ll64l for implementations in phenomenological analyses. 

Here we will use Eq. (13.431) to derive the NNLO expression c^ g and c^ q at large Q 2 ^> m 2 . 
To that end we rely on the extensions of the two-loop massive OMEs to higher orders in £ (the 
parameter of dimensional regularization) [HT1 . on the two-loop gluonic OMEs lfT2l and, most 
importantly, on the even-integer Mellin moments of the three-loop heavy-quark OMEs [|T3l along 
with their complete rcy-dependence [TT4l . The other important input consists of the results for the 
anomalous dimensions H32U33H and coefficient functions for DIS with massless quarks H17H651I661 
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Figure 2: Representative three-loop diagrams of the flavour topologies /7o2 and fl\\ (left) and //| and 
/Zfj (right) for photon-quark (left) and photon-gluon scattering (right). 

up to three loops. The necessary convolutions, Mellin transforms and inverse Mellin transforms 
are all performed with algorithms for harmonic sums WT\ and harmonic polylogarithms [1521 as 
implemented in the symbolic manipulation program FORM Il68ll . 

The coefficient functions 1 of the massless theory display one interesting feature which 
requires a comment in the context of the factorization in Eq. (13.431) . Starting from three-loop order 
C^ ht contain new quark flavour topologies, which are denoted flu and flf v with both photons 
in a squared matrix element coupling to distinct quark lines, see Fig. [21 They extend beyond the 
standard electromagnetic coupling of the heavy-quark, proportional to its squared charge e% in 
Eq. (|2.31) . as implied by the topologies //02 and //|. These new flavour classes fl\\ and fl\ x 
(being proportional to the colour factor d abc d a b c ) affect the factorization for Q 2 ^> m 2 in Eq. (13.431) 
in the following way. The available results for the massless coefficient functions in Ref. IfTTl 
(derived via the optical theorem) contain the sum of the contributions from all final state cuts 
which contribute to the inclusive massless structure function, say F 2 . However, these are not 
identical to the contributions for the heavy-quark structure function in Eq. (12.31) . the latter being 
semi-inclusive and requiring a heavy-quark pair in the final state. This is because some admissible 
final state cuts go through the quark line while others do not, e.g., the cut of the two gluon lines 
in the right diagram of Fig. |2] Currently available information is, unfortunately, not sufficient to 
disentangle the different terms. For that reason, we have decided to omit the (usually small) three- 
loop contributions proportional to flu and flf, in c} lght in the following, cf. also Appendix [Bj 

The general expansion of the heavy-quark OMEs in Eq. (13.431 ) in powers of 0C S reads 

oo co £ 

Ay = 8 ( , + £ = 8y + £ 4 £ Z>f > , (3.44) 

l=\ 1=1 k=0 

where, at each order, the terms proportional to powers of are determined by lower order OMEs 
and splitting functions (anomalous dimensions). The genuinely new £-th order information resides 
in the expressions for afj . The available results for those have been given in the literature for the 
heavy-quark mass in the on-shell scheme (pole mass) as well as in the MS scheme, see Ref. IfTBl . 
Note, that in the latter case, certain constants such as and ^2 m 2 are absent in Mellin space, 
similar to DIS structure functions with massless quarks [fPTll . 

(2) 

For the NNLO gluon coefficient function cX ' we specifically need the heavy-quark OME 

(2) 

A qh g (= AQg in the notation of Ref. lfT3l ); and for the NNLO light-quark coefficient function c 2 ' 
the pure-singlet heavy-quark OME A v q \ (= A P q in the notation of Ref. lfT3l ). The three-loop 
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expressions A^j and Ag^' ps are given in Eq. (4.39) and Eq. (4.26) of Ref. lfT3~l . respectively, with 
the //-independent terms denoted by aj^j and <3g]' pS - They can be decomposed in powers of ny as 

a Qg ~ a Qg + n f a Qg ' ( J - 43j 
(3),ps _ (3)0 (3)1 

~~ a Qq,ps +n f a Qq, P Si VAO) 



where the ny terms are now known exactly from Ref. [1141 . For the ny = terms, on the other hand, 
only a number of integer Mellin moments have been computed so far lTT3l . These are thus the only 
quantities missing for the construction of A 2 } and ci 2 j via Eq. (13.431) . 

Therefore, leaving the functions Oq®{x) and OQ g ° ps (x) unspecified for the moment, we have 
constructed the explicit results for the heavy-quark coefficient functions to NNLO. The lengthy 
exact expressions are provided in Appendix |B] see Eqs. (IB 141 ) - (1BU0I) . There are a number of 
checks on these results. First of all, up to NLO we agree with the previous results [|9l[T0l. Next, in 
the convolution Eq. (13.431) for cl 2 ] and cl 2 ] all cubic powers L? from the individual renormalization 
of the OMEs AQ g and Ag cancel, as they have to, against those from the massless coefficient 

functions C^ 1 " and C^f 11 at three loops. Moreover, all remaining powers of in and 
agree numerically with the exact results for the /j-dependence derived by renormalization-group 
methods in Refs. [|71[T5]| which are valid for all ratios of Q 2 /m 2 , if the latter are evaluated in the 
present limit Q 2 3> in 2 , see also Section^] 

As a last step, in order to arrive at phenomenologically useful results, it thus remains to use the 
available information on the Mellin moments of aQ g °(x) and Oq \ % {x) for constructing (hopefully 
sufficiently accurate) approximate expressions together with estimates of their residual uncertainty. 
This proceeds along the lines of, e.g., Refs. [[691 - 47X11 . where Mellin moments of three-loop split- 
ting and coefficient functions were successfully employed to derive useful approximation of these 
quantities prior to their exact computations [117113211331172117311 . 

A vital constraint on the functions OQ g Q [x) and Oq ° ps (x) is provided by the leading-logarithmic 
small-x behaviour of c^l and c^' determined above, see Eq. (13.411 ) and the last paragraph of 
Section [3T21 Since all other contributions to c\ ' are known, that result can be used to deduce 

L lg 

(3)0/ x _ -1, ^2 f 41984 . 160 r 224 



a Qg W = x lnx C A \~7AT + "9~ ^2 - -9- Q3 J + 0(x ) . (3.47) 
Taking into account also the double-logarithmic large- x terms, the OME a^' g ° can be written as 



fl 2°to = lA m ln ra (l-x)+/(x) + -+ (4^4+16052-22453)— (3.48) 

m=l X \ / X 

in QCD (C A = 3), where the 'smooth' function f(x) is a complicated combination of harmonic 
polylogarithms which approaches a constant for x—>l and includes powers of lnx at small x. For 
approximations based of the five known Mellin moments, N = 2, 4, . . . , 10, of Ref. lfT3l we select 
two of the large-* constants A m , the subleading small-x parameter B, and a two-parameter smooth 
function build from In 1 ' 2 * and various linear and quadratic functions in x. 
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The comparison of the Mellin transform of Eq. (13.481 ) with the known moments results in a 
system of linear equations for the chosen coefficients A m , B and the parameters of /fx). Varying 
the selected two A m and the two-parameter form of /fx), we arrive at a large number of approx- 
imations which indicate the remaining uncertainty of Oq^{x). As in Refs. 11691 - 1711 . we discard a 
small number of functions for which the system of equations is almost singular, resulting in huge 
numerical coefficients and unrealistically large oscillations of the function. In this manner we have 
determined 50 to 100 acceptable approximations and finally selected two representatives, which 
reflect the error band for most of the x-range (and particularly at small x), for the final estimate of 
agj and its residual uncertainty. 

The results of this process are shown in the left part of Fig. |3j The chosen two approximations 
(solid lines) are given by 

354.1002 ln 3 (l -x) + 479.3838 ln 2 (l -x) - 7856.784(2 -x) 
- 6233.530 ln 2 x + 9416.621x~ 1 + 1548.891 x" 1 lnx , (3.49) 
-2658.323 ln 2 (l-x) - 7449.948 ln(l-x) - 7460.002 (2 -x) 
+ 3178.819 ln 2 x + 4710.725 x" 1 + 1548.891 x 1 lnx, (3.50) 

where the coefficient of x _1 lnx is a truncation of the exact result (13 .471 ) . The average of the two 
extremes can be used as the central result of our approximation procedure. 

The same procedure has been applied to the function (x) in Eq. (13.461) . Here the lowest 
six even-integer moments have been computed in Ref. IfTBl . Taking into account that this pure- 
singlet quantity vanishes for x — > 1, our approximations are built from a subset of the functions 

(l-x)ln{ 3 ' 2 ' 1 >(l-x), (l-x)x{- 1 '°' 1 ' 2 >, xln^x, ln^x (3.51) 

together with the small-x limit given by C F /C A a^, of Eq. (13.471) . The results are shown in the 
right part of Fig. |3] and the extremal representatives have been chosen as 

(1 -x) {232.9555 ln 3 (l -x) + 1309.528 ln 2 (l -x) - 31729.716x 2 (3.52) 
+ 66638. 193x + 2825. 641X" 1 } + 41850.518xlnx + 688.396 x" 1 lnx , 

(1 -x) {126.3546 ln 2 (l-x) + 353.8539 ln(l-x) + 6787. 608x (3.53) 
+ 3780.192X" 1 } + 8571.165xlnx - 2346.893 ln 2 x + 688.396x 1 lnx , 

where, again, the coefficient of x~ ! lnx arises from the truncated exact result, and the average of 
Eqs. (13.521) and (13.531) provides the central result. 

Using these approximations, the high-scale coefficient functions c 2 g and c 2 q can now be as- 
sembled up to NNLO using Eqs. (1B14I) - (1BU0I) . The successive LO, NLO and NNLO results are 
shown in Fig.@]at Q 2 = 100 GeV 2 , a scale that is low enough to be relevant to the HERA measure- 
ments and sufficiently high for the safe applicability of the Q 2 ^> m approximation (13.43b . We 
observe a satisfactory convergence, especially of the dominant quantity c 2 g , over a wide range of x, 
together with the well-understood large-x region of soft-gluon enhancement for c 2 , cf. Ref. f|42]]. 



a {3)0 (x) 
a QgA x > 



a (3)0 fx) 



a 



(3)0 
Qq,ps,B 



to 
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Figure 3: A large number of approximations for the OMEs (left) and a>Q q ps (right) compatible with 
the leading-logarithmic small- x behaviour (13.471 ) and the known first five even-integer Mellin moments ifTBl . 
The solid curves show the respective selected representatives ( 13.491 ), ( 13.501 ) and ( 13.521 ), ( 13.531 ). 
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X X 

Figure 4: The perturbative expansion of the coefficient functions c 2 (left) and c 2 (right) for heavy-quark 
DIS at a large scale Q 2 S> m 2 , with the respective uncertainty bands due to the OMEs in Eqs. ( 13.491 ), (13.501) 
and (13.52b . (13.53b . The results are shown for three light flavours and a typical value of a s at fi 2 = Q 2 + 4m 2 . 
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For both coefficient functions the uncertainty band due to the heavy-quark OMEs Aq s and Aq 
as given by Eqs. (13.491 ), (13.501) and (13.521) . (13.531) . respectively, is phenomenologically relevant only 
at small x, up to x ~ 10~ 2 . The uncertainty in this region is due to that of the (large) coefficient 
of the 1 jx terms in Eq. (13.481) and its pure-singlet counterpart. The computation of more Mellin 
moments along lines of Ref. lfT3l would certainly help to improve the approximation also here; 
however the issue would be settled for all practical purposes by extending the results of Ref. [[HI to 
the next-to-leading small-* terms at order a s 3 . 




Figure 5: Top panels: the coefficients of the leading small -x / high-q logarithm for the contributions c 2 g , 
£ = 0, 1,2, to the NNLO gluon coefficient function defined in Eq. (I2.6I ). Middle and bottom panels: the re- 
spective next-to-leading q° coefficients for and c^J . The solid (black) lines are the exact all- Q 2 results, 
the dashed (blue) ones the high-scale asymptotic results; the dotted (red) low low-scale extrapolations are 
discussed in Section 0] below. Also illustrated, by the thin (green) lines in the bottom panels, is the small 

(2) (2) 

next-to-leading high-q deviation of cX from the 'Casimir-scaled' results for c) . 
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A crucial question is, obviously, down to which values of ^ = Q 2 /m 2 these high-scale asymp- 
totic results are applicable. For the important high-energy (small-*) contributions, which one may 
expect to be least affected by m 2 / Q 2 corrections, this issue is addressed in Fig. |5]for the quantities 
, i = q,g, in Eq. (12.61) . Here the exact leading, cf. Eqs. (13.391 ) and (13.411) . and next-to-leading 
high-ri results are compared to the respective high-scale expressions. Except for the L 2 quantities 

(2 2) 

c\ ■ in Eq. (12.61) fixed by LO information, the asymptotic results provide good approximations of 

the lnr| parts of c^f^ and the lnr| and T|° parts of c?! down to £, ~ 4, but then deteriorate dra- 
matically towards % = 1 . This behaviour can be used to extrapolate the (unfortunately still rather 
uncertain) asymptotic coefficients of r|° for c^j to small \. This extrapolation will be discussed 
in the next section, but the corresponding curves have already been included in the figure. 



4 Approximate coefficient functions at NNLO 

We are now in a position to construct improved NNLO approximations for the genuinely new parts 
of the NNLO gluon and pure-singlet (ef t ) coefficient functions, c^g and , for heavy-quark 
DIS. Specifically, we will discuss how to combine the above (approximate) expressions in the three 
kinematic limits, i.e., threshold, high-energy and large Q 2 , in order to arrive at functional forms 
which smoothly interpolate over the full relevant range in T| and £ (recall, again, Eq. (12.4b ). 

All logarithmic terms proportional to powers of L„ in Eq. (12.61) are known exactly from stan- 
dard renormalization- group methods. For completeness, we give their explicit form in the MS 
scheme. Suppressing all arguments, these coefficients read up to NNLO H7l [T5l[2TI 
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where, as in Eq. (13.431 ) above, (g> denotes the standard convolution. The splitting functions pfj' (x) 
can be taken from [|32~ll33l and (3o and (3i are the standard coefficients of the QCD beta-function, 
normalized such that |3o = 1 1/3 Ca — 2/3 «y. The powers of 471 account for the normalization 
of Refs. [@1[51 adopted in this article, cf. Eq. (12.61) . Lacking more exact third-order information, 
in particular c 2g in Eq. (14.51) will provide important guidance below for the assembly of the 



(2,0) 



approximate NNLO result for the most important quantity c 2 g 
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Two issues need to be considered for the combination of the results from the various kinematic 
regions. Firstly, one has to merge the available information from threshold at low P (and r|) 
with the high-energy limit at large T|. Since we use the complete Born coefficient function (12.71) 
which vanishes for r\ — > °° as the prefactor in the threshold expansion (13 .18b . this can be done by 
multiplying the high-energy terms by sufficiently (but not excessively) high powers of |3. Secondly, 
the asymptotic expressions of the large-S, limit need to be joined with the low-£, region. To that end 
we employ the function 

which provides a smooth transition between these two regimes. The parameters 2 and 4 in the 
exponent are chosen in such a way that the transition is sufficiently rapid at values close to ^ = 5. 
This choice is motivated by the above discussion of Fig. [5] and the finding of Ref. [91 that, e.g. for 
C2g \ the asymptotic limit (IBI6I) represents the exact result with a high accuracy already at ^ > 10. 

These considerations lead the following ansatz for c^g (cf. also the NLO parametrizations in 
Ref. (51) which we test below, 

C £ - 4r + (w(Q) pM!i asm + nw 3 K )LLx cw) • (4 - 8) 

Here c^g an ^ c^ 1 ^™ denote the threshold expansion in Eq. (13.161 ) and the asymptotic result 

in Eq. (IBI6|) . respectively. The leading term in the small-x limit in Eq. (13.381) defines c^ 1 ^* and 
provides the high-T] tail in the low-£, region which is smoothly matched with the factor rp/ (C+r^) . 
The values of y, C and k in Eq. (14.81) will be specified below. 

(2) 

Correspondingly, for c\ g we use 

cfl K C PJ* + (l- /S ))p*4^» +/S) p3(_4^xg + c (2J^L^ (49) 

with c^g thl and c^ 3 ™ of Eqs. (13.181) and (1B18I) . The description of the low-^ region in Eq. (14.91 ) 

deserves particular attention. Here c^ LLx is the leading contribution (13.391 ) in the small-* limit, 
and we have divided out the factor \nx in order to be able to substitute lnx — > — hit], which is 
valid at high energies and determines the slope in r\ at all values of ^. The next-to-leading large-T] 
term, denoted by c 2 2 ] NLL in Eq. (14.91 ), is related (but, due to r\ = t,/ {Ax) + 0(1) not identical) to 
the 0(x°) contribution in Eq. (I3.39I ). It is currently unknown in the low-^ region. We will derive 

constraints on c 2 2 ^ NLL below, although, even at large ^, we still have to cope with the uncertainty 
of the heavy-quark OME Ag g estimated by Eqs. (13.491 ) and (13. 501 ). Recall also that the unknown 
next term in the low-T] threshold expansion (13.181 ) consists of terms proportional to the Born result 
as discussed in Section I3T1 

Therefore we adopt the following strategy for applying Eqs. (14.81) and (I4.9I ). When merging 
the contributions from the various regions, we account for the residual uncertainties and construct 
two enveloping curves designed to span the uncertainty band in the entire kinematic plane. Thus 
at NLO we define 
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c 0,0M = c M) fa _ c (Wco rat+(1 _ /K))pc a0)«, m+/fi)p 3^)LL,_^^ (4]0) 
C »K0,, S = c ,,.0„ to +(1 _ /ffi)) p3 c (l 8 0).m +/ffi)p 3^0)LLx_ni 5 ^ (4H) 

Here c^ 1 ^"" and C2 ^°) asm are given by the L M -independent parts of Eqs. (13.161) and (IBI5I) . while 

the term c ^°) const subtracted in Eq. (14.101) summarizes the NLO constants in (3 in Eq. (13.171) . The 
suppression parameters take the values 

y=1.0, C = 42.5, and 8 = 0.8, D=19.4. (4.12) 

They are needed in order to suppress the coefficients of r\° at low r| and have been determined by 
fitting the high-r| tail of c^g using the MlNUlT package of the CERN FORTRAN library. It is not 
possible to achieve agreement with the exact result without such an additional suppression factor. 

The resulting two approximations are compared, after adding the exact scale-dependent con- 
tribution, in the left parts of Fig. |6]to the parametrized [0 exact result of Ref. flU at four typical 
values of £, for the standard choice /j 2 = Q 2 + 4m 2 of the (identical) renormalization and factor- 
ization scale. It is worthwhile to note not only that the approximations (14.101) and (14.1 II) indeed 
span the exact results, but also that there is, at this scale, a considerable cancellation between the 
high-energy constants in c^ 1 ^ and ' '• 



At NNLO we proceed in a similar manner. As indicated below Eq. (14.61) . we first consider the 

(2 1) 

known quantity c\ ? and define 



(2,l),A _ (2,1) thr n wtun (2,l)asm frn R 3 f J2,1)LLx lnr| (2,1)NLL if 
C 2,g ~ C 2,g + l 1 MSJJP C 2 ,g + J c 2,g \nx +C 2,g C+Tl 

4f' B = 4f thr + (i-^))P 3 4! ? 1)asm +^)P 3 (-4f LLx S+^ 



(4.13) 

(2,1) NLL r| S 
8 D+if, 

(4.14) 



with c^^ thl , cf^ g l ^ asm and c^^ LLx respectively given by terms linear in in Eqs. (13.181) . (1B18I) 

and (13.391) . For the next-to-leading term at small-* denoted by c^' g l ^ NLL we can derive from 
Eq. (14.51) the parametrizations 

c 2 2 g 1)NLL (^) - 0-16086 - 0.00711 ln^ - 0.00549 ln 2 ^ . (4.15) 
The damping of this NLL high-energy tail can be approximated by (again using MlNUlT) 

y=1.0, C = 20.0, and 8 = 0.8, D=10.7. (4.16) 

These approximations for ^ are compared in the right parts of Fig. |6]with the exact result (14.51) . 
for which these graphs supersede previous studies in Refs. [0[3O]|. Also these approximations are 
completely satisfactory, taking into account that a term corresponding to c ^°) const m gq (|4. 1Q|) 
has not been included in Eqs. (14.131) and (14.141) . 



22 



0.6 

0.4 

0.2 



0.4 
0.3 
0.2 
0.1 


-0.1 

0.2 
0.1 


-0.1 



llllllllj llllllllj I lllllllj llllllllj llllllllj llllllllj llllllllj lllllllll 



t> = 2 




i i 



J I 



i i 



J I 



i 1 1 



_ 1 II lllllj 


1 1 1 lllllj 


i iiiinij i ii iiiiij i iiiinij 


lllllllj 


1 1 1 lllllj 


1 1 nil 








I 


= 5 






J? 


\\ 












\\ 






Ml jlll 






w 








1 1 1 






— 




i i i 


~ 1 


1 


1 1 1 


1 


1 


1 1 nil 



0.15 
0.1 
0.05 


-0.05 
0.05 



- 



J lllllllj llllllllj llllllllj llllllllj llllllllj llllllllj llllllllj I MINI 

c 2 (2 g 1} (ti4) \ = 2 




1 1 1 lllllj 


1 lllllllj 1 lllllllj 


1 lllllllj 1 II Mill 


1 lllllllj 1 lllllllj 1 


II Mill 








§ = 10 






/£' 


*\\ 










\ 










\ 












\t^y^ 




1 


1 1 


1 


1 1 I 


II Mill 



-0.05 
0.06 

0.04 

0.02 



-0.02 
-0.04 



1 lllllllj 1 lllllllj 1 lllllllj 1 lllllllj 1 lllllllj 


i iiiinij i iiiinij i iimn 










exact 


\ * 

\\ - 

\ • - 

\\ 


: approx AB 

i i i i i 


\ '- 

\\ 

1 1 n 



1 1 1 lllllj 


1 lllllllj 


1 1 1 lllllj 1 1 1 lllllj 


lllllllj 


1 1 nun 


i iiiinij i 


lllllll 






*\ 




t 


= 10 








\\ 
\\ 
























1 


1 


1 1 


1 





1 1 


iimn 



0.2 - 



0.1 - 







llllllllj llllllllj llllllllj llllllllj llllllllj llllllllj llllllllj I lllllll 




0.02 - 



I I I I I I I 



1 1 1 lllllj 


1 lllllll 


1 1 1 lllllj 1 1 1 lllllj 


1 lllllllj 


1 1 1 lllllj 


1 lllllllj 1 


iiniii 










\ 


= 20 




































N \ 














v\ 
























1 





1 1 


1 


1 


i i 


Miiiii 



10 



10 



10' 



10 

,(2,1) 



10 



10' 



Figure 6: Comparison of the exact result for c\ l, g and c^' g l> , at four representative values of t, = Q 2 /m 2 , 
with the approximations based on the threshold, high-r) and high-<2 2 limits as respectively specified in 
Eqs. (14.101 ) and (14.111) and Eqs. (14.131) and (14.141) . The first-order results are shown at a physically relevant 
scale jj after adding the exact scale-dependent contribution to all three curves. 
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Encouraged by these results, we now apply the same procedure to the main new NNLO coeffi- 
cient function ci 2 ' '* and write 

(2,0),A _ (2.0) thr n ff p\\a „ (2,0)asm,A 
c 2,g - c 2,g + U ASJJP C 2ig 

+ /© p 3 ( - ,g 0,LLx S + ^ 0)NLM c^) <4.n) 



and 



cg 0) ' S = 4 2 / )thr + /(^)24y )con V(l-m))(3 3 cg 0)asm ' S 

+ /(SJP ^ _c 2,g i„ jc + c 2, g D+r ,8 J (4- 1 8) 



with c^ g °^ tbl and c ^°) asm given by the -independent terms of Eqs. (13.181 ) and (IBI8I) . It is un- 
derstood that in the latter case the results for tfgg A an( * a gg e m Eqs. (13 .49b and (13.501 ) are used to 

account for the uncertainty band due to the heavy-quark OME Aq s . The term c ^°) const denotes 
the scale-independent part of the Pade estimate Eq. (13.191 ) of the NNLO constant in (3. The factor 
of two in front originates from assigning a 100% uncertainty to this estimate, i.e., we add zero in 
Eq. (14.171) and 2c2 2 g °^ const in Eq. (14.181) . The values for the suppression parameters y, 8, C and D 

are taken from Eq. (14.161 ). and the functions cj 2 ^ NLL in the low-^ region at high-r] are estimated 
via the error band 

c (2,0)NLL,A ( ^ = Q ^lnjy _ Q 2g ^ (4 J9) 

cJ 2 / )NLM (^) = 0.055 (j^Y - 0.423 . (4.20) 

These are the extrapolations already shown in Fig. [5] They smoothly continue the large-scale 
results at 5 < \ < 10 to the lower value of \ considered here, \ = 1, in a manner suggested by the 
corresponding low-^ behaviour of and c^g ■ The fine details of Eqs. (14.191) and (14.201) do not 
matter, of course, due to the large uncertainty of the underlying high-S, results (13.491) and (13.501) . 

In Fig. [7] we display these approximate NNLO results for c^g for seven values of £, in the 

(2 1) (1 21 

range 1 < £, < 100. The known scale-dependent terms with c\ ^ and c\ g are then added in 
Fig.[8]to these expressions for the standard scale /j 2 = Q 2 + 4m 2 which we will also use in Section[5] 
below. The lack of a stronger constraint on the next-to-leading high-T] coefficient leads to a large 
uncertainty in particular at low where it extends down to T] « 10 in the construction detailed 
above. We have investigated other approaches to this important mid-^ region, but have been unable 
to derive a reliable stronger constraint with the available information on c^"^ . The present NNLO 
results in Fig. [8]are consistent, albeit if with a large uncertainty, with the 'high-energy censorship' 
already seen for c^l in Fig.[6l i.e., a small contribution of the high-energy tail to the total coefficient 

(2 i) 

function despite large contributions to the individual components c\ ' , £ = 0, 1, 2. 
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Figure 7: The two representative approximations (I4.17t and (I4.18l l for the 'irreducible' scale-independent 
contribution c^"p to the NNLO gluon coefficient function for the structure function F 2 in heavy-quark 
production in DIS in Eq. ( 12.21 ). shown for the experimentally most relevant range of ^ = Q 2 /m 2 , where m 
denotes the heavy-quark mass, for = 3 light flavours. Also shown, for £, < 5, is the low-^ improvement 
of the approximation 'B' as discussed at the end of this section. 



25 



iiiiiiiil iiiimij i iiiiiiij iiiiinij iiiimij iiiiiiiij iiiiiiiii iiiiimi 

5=1 




iiiimij iiiiinij iiiiiiiij iiiimij iiiiinij iiiimij iiiiiiiij iiiiiiiii 

t> = 2 




_Q 2 I I I I I I I I I _Q 2 1 L 

iiiiiiiii iiiimij i iiiiiiij iiiiinij iiiimij iiiiinij iiiimij iiiiiij iiiiiiiiij iiiiinij iiiiiiiij iiiimij iiiiinij iiiimij iiniiiij iinn 




5= 10 



0.3 
0.2 
0.1 


-0.1 F~i I I I I I I I i niiiiil -0.1 Pi I I I I I I I i i liTiii 

0.2 |_ IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIJ 0.2 [jm™ Illlllllj IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII IIIIIIIII I Mill 





lllj I Illllllj 



5 = 50 -_ 




q Q5 r I I i i i i i in 

-4 -2 2 

10 10 1 10 

cfW), |i 2 = Q 2 + 4m 2 



ml I 



0.06 
0.04 
0.02 


-0.02 



10 



1 1 1 II lllj 


1 Illllllj 1 II Mill 


1 Illllllj 1 Illllllj 


1 Illllllj 1 Illllllj 1 II Mill 








5 = ioo -_ 














\ 








\ 












~ 1 


1 " 


1 "I 


1 1 iii 


)- 4 


io - 2 


1 





approx. AB 
approx. B , {low ^} 



(2) 

Figure 8: As Fig.|7J but for the complete coefficient function c\g, obtained by adding the known scale- 
dependent parts, for a standard choice of the renormalization and mass-factorization scale n 1 = fij = ^. 
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(4.21) 



For completeness, we apply the same procedure to the pure-singlet coefficient function c 2 q 
which exhibits a closely related high-energy behaviour, recall the end of Section lX2l and the bottom 
parts of Fig. [5} but no low-T) threshold enhancement. Suppressing the quantities corresponding to 
Fig.[6]for brevity, we directly turn to the new NNLO quantity c^ q °^ which thus we approximate by 

(2,0),A _ , f(%\\a J2,0)asm,A ff *s R 3( (2,0) LLx lnT| (2,0)NLL,A 

C 2,q ~ U -/ISJJP C 2,q [~ C A C 2,g Inx + C 2,q C+T|Y 

(2,0), B n ff .,s r3 (2,0) asmB , f R 3 f C F (2,0)LLx lnr| (2,0)NLL,fi r| 8 

C 2,q ~ y L -JVV)P % q +/^JP Ca C 2g \ nx + c 2,q D+r ,8 

(4.22) 

The high-£, asymptotic result c ^°) asm i s given by the //-independent part of Eq. (|BI10|) . which 
presently includes the two approximations (13.521 ) and (13.531) for the pure- singlet operator matrix 
element a^j^. The values Y, 8, C and D are again taken from Eq. (I4.16I ). Analogous to Eqs. (14.191 ) 
and (14.201) . the low-^ extrapolations of the high-T] constant, constrained by the 1 /jc-coefficients in 
Eqs. (13.521) and (13.531) . are chosen as (see again also Fig. |5]) 

c (2,o)nll.a^ = o.004 (^|V - 0.125 , (4.23) 

c (^0)nll, S( ^ = Q245 /ln|\ 2 _ Q 1? (4 24) 

The resulting approximations for and the (small - compare the scales of the right parts to 

those of Fig. [8]) complete coefficient function c\ q at// 2 = Q 2 +4m 2 are illustrated in Fig. |9] 

So far we have separately considered approximations to the NNLO gluon and pure-singlet 
quark coefficient functions. However, as mention above, the small-x/large-r) limits of these two 
functions are closely related; and while their C F /C A relation does not hold exactly beyond the 
leading logarithms, it still approximately holds for the NNLO T| -terms as shown already in the 
bottom panels of Fig. [5] in Section [331 In particular, the numerical breaking of this relation for 
these terms is much smaller than the current uncertainties of c^f^ and ci 2 ' '' in the low-£ region. 

This opens up the possibility of one final improvement, i.e., using the C F /C A relation and the 
better constrained coefficient function to reduce the uncertainty of its more uncertain counterpart. 
For both the OME and the coefficient function the high-energy uncertainty is smaller in the quark 
case, as can be seen from Figs.[3]and[5] which is unsurprising given that one more high-^ moment 

is known for the OME in this case lfT3l . which is moreover an 'easier' function with a^} Q „(x) — > 
for x —7- 1. Since the 'A' approximations are very close to the C F /C A relation, cf. Fig. [5] this leads 
to replacing Eq. (14.201) by 

c (2,0)NLL,fi' = £4 c (2,0)NLL,S ^ ^ 

2^8 Cf 2,q 

with the right-hand side given by Eq. (14.241) . The resulting improved low-^ estimate has been 
included in the first three bins in Figs. [7] and [8] and it is this reduced uncertainty estimate that we 
will use when we study the approximate NNLO results for F 2 charm and i^ bottom in the next section. 
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Figure 9: As Figs. [7] and [8j but for the pure-singlet (light-) quark coefficient function c^q in terms of the 
approximations (14.211 ) and (14.221 ) at some low and intermediate values of \. 
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5 Charm and bottom structure functions at NNLO 



We finally illustrate the impact of the above approximate NNLO coefficient functions on the heavy- 
quark contribution to the structure function F 2 in Eqs. (12.21) and (12.31) . Since present-day parton 
distributions are rather well constrained for most of the kinematic range relevant to the charm and 
bottom structure functions, we confine ourselves to one set, that of Ref. 11741 . See Refs. H63U75I - I791 
for other recent PDF fits including NNLO corrections. For the heavy-quark pole masses we use 

m 2 = 2GeV 2 and m b = 4.5 GeV , (5.1) 

and our default choice for the renormalization and mass -factorization scales is 

li 2 = fij = = Q 2 + 4m 2 h , h = c,b. (5.2) 

In Fig. [10] the resulting charm structure function F 2 is shown for seven scales Q 2 which cor- 
respond to the ^-values chosen in Figs. 7-9 above. The dotted and short-dashed curves are the 
contributions of the LO and NLO coefficient functions to the NNLO results, i.e., these quantities 
have been calculated using the NNLO values of the strong coupling and the parton distributions 
as determined in Ref. ir74~l . The solid and dash-dotted curves show the NNLO corrections for the 
approximations 'A' and 'B' constructed in the previous section, including the improvement (14.251 ) 
of the latter. Where the approximations are sufficiently accurate, roughly speaking at x > 10~ 3 , the 
overall effect of the NNLO coefficient functions is positive but considerably smaller than that of 
the NLO contribution, indicating a good convergence of the perturbation series. At x < 1(T 3 the 
present results are inconclusive in this respect in particular, unfortunately, in the important low- 
and medium-<2 2 region. Also shown, by the long-dashed curves, are the physical NLO results 
involving the NLO coupling constants and parton distributions according to Ref. [1741 . 

The corresponding results for the bottom-production contribution to F 2 are illustrated in Fig.fTTI 
for the slightly shifted region 0.5 < ^ < 50. The pattern of the corrections is qualitatively similar 
to that for charm production in Fig. [10] Again the NNLO approximation uncertainty is largest at 
x < 10~ 3 and 1 < £ < 5, here corresponding to 20 GeV 2 < Q 2 < 100 GeV 2 , where a pure threshold 
approximation can be insufficient and the mid-r| constraints of the high-scale expressions are not 
applicable yet. Note that this figure includes a bin at Q 2 < m\, where the impact of the medium- 
and large-T] parts of the coefficient functions is smaller than at higher scales. In general, the size 
and uncertainty of the NNLO corrections are smaller for F 2 than for F 2 at the same values of £ 
due to the smaller values of oc s and the steeper small-* PDFs at the corresponding higher scales. 

Returning to the case of charm production which has a far larger impact on the determination 
of especially the gluon distribution from DIS data, Fi g s . ITU and [T3l provide a more detailed look at 
the NNLO results and their remaining uncertainties for two typical values of Q 2 . The former figure 
shows the dependence on the scale /u for the range m c < /u < 2/j st i with ^ st( j = /u of Eq. (15.21) which 
includes both hard scales entering the problem, 2m c and Q. Except at very low scales, ^ < 2m c , 
(and possibly at very small x) the NNLO scale variation of F 2 is smaller than its NLO counterpart 
given at /j = ^ st d by the long-dashed curves in Fig. [T0l This is also obvious from Fig. [T21 where the 
NLO and average (j( A'+'B')) NNLO scale uncertainties, estimated via the maximal and minimal 
values in the conventional interval [j /J s td ? 2^ st d], are compared with the present approximation 
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Figure 10: The charm-production structure function F 2 c (x,<2 2 ) in the experimentally important region 
2 GeV 2 < Q 2 < 200 GeV 2 . The effects of the LO, NLO and approximate NNLO coefficient functions are 
illustrated by using in Eq. (12.31 ) the NNLO values of a s and the parton distributions of the recent ABM 1 1 
fit 11741 . Also shown are the NLO results obtained for the corresponding NLO values of a s and the PDFs. 
All curves are shown for a pole mass of m 2 = 2 GeV 2 and the standard scale \i 2 = Q 2 + 4m 2 . 
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uncertainty which becomes important at x ~ 10~ 3 . The relatively small NNLO improvement at 
large x indicates the need for yet higher orders. Here a first step would be to include the (almost 
complete) NNLL threshold resummation discussed in Section |3T| 

6 Summary and outlook 

The production of heavy quarks, in particular charm, contributes a sizeable fraction of the proton 
structure function F 2 measured at HERA, thus affecting fits of the parton distributions and hence 
the predictions for hard processes at other colliders such as the LHC. Unlike the case of massless 
quarks, the corresponding partonic cross sections (coefficient functions) are not fully known yet at 
the next-to-next-to-leading order (NNLO) of perturbative QCD. In this article, we have collected 
and extended the partial NNLO results for three kinematic limits, i.e., the threshold region with 
s > 4m , the limit of high partonic CM energies, s ^> m 2 , and the high-scale region Q 2 3> m 2 , and 
then combined these results into approximate expressions covering most of the kinematic plane. 

(2) 

For the NNLO gluon coefficient function c 2 g we have determined all threshold contributions 
which are logarithmically enhanced in the heavy-quark velocity p\ the so-called Sudakov loga- 
rithms, as well as the Coulomb corrections leading to powers of 1 / (3 (Eq. (13.181) in Section. [37X1) . 
In the high-energy limit, we have employed results of the corresponding ' small- x' resumma- 
tion [8] to derive explicit expression for the (closely related) leading-logarithmic In s contribution 
to the gluon and quark coefficient functions at order a 3 (Eqs. (13.391) and (13.411) in Section I3T21) . 
Finally, in the high-scale limit where m 2 / Q 2 power corrections can be disregarded, we have utilized 
the mass-factorization formula to derive exact expressions for the NNLO gluon and (pure- singlet) 
quark coefficient functions c 2 ^ and c^ q (Section [331 and Eqs. (1BI7I) and (1BU0I) in Appendix B). 

The latter results are complete for the contributions proportional to ny, the number of light 
flavours, and rely on the three-loop operator-matrix elements (OMEs) of Ref. lfT4ll for the Q 2 
independent contribution. Their counterparts for the /^-independent contributions are not fully 
known yet, though, so we had to rely on the low even-iV Mellin moments computed in Ref. fP3l 
and our above small-* result for deriving the approximate expressions given in Eqs. (13.491) - (13.531) 
in Section 13.31 The largest uncertainty in these results is due to the next-to-leading small-* (at 
NNLO behaving as 5° at large s) contributions which are not tightly constrained at this point. 

We have then combined the results of the various regions to provide the best possible approxi- 
mations for the coefficient functions for heavy-quark production in DIS at all values of s and Q 2 , 
carefully modeling and estimating the impact of missing information such as the NNLO threshold 
constant and of the uncertainty of the high-5 constants and their dampening towards medium val- 
ues of s. The corresponding results are given in Eq. (|4.171 i - (I4.251 i in Section HI The implications 
of these uncertainties for the charm and bottom structure functions F 2 c (jc, Q 2 ) and F 2 (x, Q 2 ) have 
been illustrated in Section|5l where we have also addressed the dependence on the renormalization 
and factorization scale /j which can signal the importance of yet higher orders. The uncertainty 
due to the missing information on the NNLO coefficient functions is large at x < 10~ 3 especially 
at low and medium scales Q 2 , but irrelevant at much larger values of x. The dominant uncertainty 
in that region (where the absolute values of F 2 C and F 2 are much smaller) is due to higher orders, 
as shown by the rather modest NNLO improvement of the /j-dependence. 
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At high scales, Q 2 > 10 m 2 , the remaining small-* problem can be solved by extending the 
fc r -factorization results of Ref. [O to the next-to-leading high-energy contributions at order a s 3 or 
by extending the calculations of |fl4]| to the n^-independent OMEs. Already the extension of the 
fixed Mellin-iV computations of Ref. |[T3l by a couple of moments, however, would facilitate a 
considerable reduction of the uncertainties. These calculations could also be employed to con- 
siderably reduce, but presumably not remove, the large small-x uncertainties at low and medium 
scales, Q 2 < 10m 2 . It is also worthwhile to note that there is a close relation between the next- 
to-leading high-energy contributions in the present case and in heavy-quark hadro-production, see 
Refs. [|8l[50l. This has been exploited in Ref. [80J to estimate the high-energy tail of the coefficient 
functions for top-pair production at the LHC, which complements numerically determined results 
on NNLO heavy-quark production in hadron-hadron collisions lf8TI valid in the region where the 
partonic CM energy s is not too large, s/4m 2 < 100. 

The uncertainties at low Q 2 may be removed to a fully satisfactory extent only by an exact third- 
order calculation of heavy-quark production in DIS. This is a formidable task, including three- 
loop Feynman diagrams with two scales, i.e., an internal mass and an off-shell leg. A useful, but 
conceptually easier intermediate step would be a calculation in the threshold limit determining the 
missing threshold 'constant' (i.e., the contribution proportional to the Born cross section) at order 
a 3 . This result would also provide the only piece missing for extending our results in Section l3Tl to 
a NNLO+NNLL (next-to-next-to-leading logarithmic) threshold resummation. This resummation, 
even without the NNLO constant, provides the natural framework for a first step beyond the third 
order, e.g., for addressing the still rather large scale uncertainty at large x. 

At this point, though, our approximate third-order results represent the best and most complete 
predictions of heavy-quark DIS in the important context of NNLO determinations of parton distri- 
butions from inclusive and heavy-quark structure functions measured at HERA. Given the wider 
implications of those analyses on LHC physics, however, further improvements in particular at 
small x would be extremely welcome. For use until then, Form files and Fortran subroutines 
with our main results can be obtained from the preprint server fhttp: //arXiv. org| by download- 
ing the source of this article. Furthermore they are available from the authors upon request. 
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A Gluon coefficient function at NLO for s — > Am 2 



In this first appendix we derive the functions c (£) and c (^) in Eqs. (13.101) and (13.1 II) for the 
gluon coefficient function c 9 „ at NLO. This derivation starts from the coefficient function in the 
differential kinematics of one-particle inclusive DIS D4J, 



£(1+P) f s^ d 2 4 n l( S ',t h U h l? 

d(-ti) / ds 4 — ,g 
t(i-p) 7o 



<i?i ds4 



(EH) 



where the kinematic variables for the subprocess + J*(q) —> 1h(pi) +X[lh] are defined as 
*l = (zp ~~ Pi) 2 — m2 ' M l = (q — Pi) 2 — m 2 and s' = s + Q 2 . The variable s 4 = s' + t\ + u\ = 
M\ - m 2 > measures the inelasticity of the subprocess, with s 4 — > corresponding to the limit 
of soft gluon emission. The upper limit of s 4 is given by 
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The LO differential coefficient function is given by [|4l[T8l 
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The NLO differential coefficient function was calculated in Ref. H. In the threshold region s 4 ~ 0, 
it can be simplified as 01EIESD 
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and the +-distribution in 54 defined as 



'hf(s 4 /m 2 Y 
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The conversion of to the standard MS scheme for a s proceeds via the well-known relation 
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and affects the function co(^) for the scale-dependent matching constant, see Eq. (IAI14I) below. 
The scale-independent function R(s', ti,u\) in Eq. (IA[6I) originates from the soft- and virtual- gluon 
contributions, for which the analytic expression is quite lengthy and available only via a FORTRAN 
code by the authors of Ref. [|4]|. 

In order to identify the leading contribution for s -+ Am 2 , i.e., in the limit [3 — > 0, we first 
consider Eq. (IAI1I) for the LO coefficient function with the integrand given by Eq. (|A13I) . After 
the trivial 54-integral, we have u\ = —s' — 1\, and the leading contribution in the fi -integral can be 
obtained by the simple replacement 



-+-2' 



and multiplying with an overall factor s'$. This can be understood by observing 
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where the first term on the right-hand side gives the leading contribution in the |3 unless the deriva- 
tive of f(t\) at t\ = —s'/2 contains negative powers of p\ Furthermore, we perform expansion in 
P such that s = 4/« 2 (l + 0((3 2 )) and we recover from Eq. (|AI41 i the threshold approximation of the 
LO coefficient function in Eq. (13.151) . 

At NLO we first note that constants and positive powers in 54, which are ignored in factorization 
ansatz of Eq. (IAI5|) . yield corrections of higher order in P after integration due to s , 4 nax ~ C((3 2 ). 
Then the replacement of (1AI101 ) can effectively be applied to a Bom in (IAI5I) . so that the 54-integrals 
become trivial with the +- distribution from Eq. (1AI81 ). Care needs to be taken in the subsequent 
^-integration, though, if we cannot rely on the replacement rule (|AI10|) . This exception only arises 
if the terms \n n (s™ ax /m 2 ) (n= 1,2) appear via Eq. (IAI8I ). because the coefficients of their Taylor 
expansion at t\ = — s'/2 contain negative power of p\ thus invalidating the assumptions underlying 
Eq. (1AI11I) . These terms give rise to the threshold logarithms, and we obtain by explicit calculation 
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Using these formulae together with Eq. (1AI101 ) and the expansions in p\ we recover the threshold 
approximation of the gluon coefficient function in Eq. (13.161) . now with explicit results for the 
the constant terms c (£) and c (^) at NLO. 

The scale-dependent term cq(4) which is also fixed by renormalization- group constraints can 
be easily read off from Eq. (IAI6I ). 



AC, 



2 + ln 1 + 



©14) 



where the last term is due to the transformation of the strong coupling from the MOM scheme in 01 
to a s in the standard MS scheme lfT2ll62ll . cf. Eq. (IAI9I) . The calculation of the scale-independent 
term co(^) on the other hand requires some automated manipulations of R(s',t\, u\), based on the 
Fortran code of flU, which finally lead to the new result (13.101 ) in Section |3H 



B Exact results at asymptotic values Q 2 m 2 

We close by presenting the exact expressions for the heavy-quark coefficient function H2 )g and 

at high scales, Q 2 3> m 2 . In order to shorten the results, we have always separated the contribution 

(£) 

of the coefficient functions with massless quarks in photon-exchange DIS, denoted as c\ I below, 
as computed in Ref. [fTTll at £ < 3. See also Refs. J65l[66l for earlier results up to the second order. 

For easier comparison with the literature, we will use the following established normalization 
of these asymptotic coefficient function, cf. Ref. llMTOll. 

H 2 , g (x,Q\m 2 ) = a s H^ g (x,Q 2 ,m 2 ) + a 2 H®(x,Q 2 ,m 2 ) + a* H®(x,Q 2 ,m 2 ) + 0(a 4 s ) ©1) 

and correspondingly for H^. As above we expand in terms of a s = oc s / (4ft ) and, again following 
Refs. [|9l[T0l. this expansion is performed in the MS scheme with a s (nr+ 1). All terms originating 
from subsequent matching a s (n^ + 1) — > OL s (n^) with the help of the decoupling formulae ll60~ll6Tl 
will be presented separately as well. The relation of, e.g., the quantity (IB 111) to the gluon coefficient 
function c 2g , defined in Eq. (12.31) in the normalization of Refs. flUEl, is given by 

c 2 , g 0U) = ^{//^(*,eW) + a s /^ dBj2) 

where we have suppressed the additional dependence on /j = /u^ = /u r . 

The explicit results are presented in terms of harmonic poly logarithms H mi m w ( x ) with m j — 

0, ± 1 , where our notation follows Ref. [|52l to which the reader is referred for a detailed discussion. 
For chains of indices zero we again employ the abbreviated notation 

H 0,...,0,±1,0,...,0,±1, ...(-*) = H ±(m+l),±(n+l),... © 3 ) 
m n 

in which also the argument of has been suppressed for brevity. The numerical evaluation of 
the the harmonic polylogarithms relies on the Fortran package of Ref. If82l and its weight-five 
extension provided by the authors. 
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The leading-order coefficient H^g in the expansion (|BIH) reads 



H { 2 %,Q 2 ,m 2 ) = 2(l-2x + 2x z )L Q + c^(x,n / + l)/(n f + 1) 



dj4) 



where c^j (x, + 1 ) denotes the one-loop gluon coefficient function with massless quarks at /j = Q 
as given in Eq. (B.4) of Ref. lTT7ll . and Lq had been defined in Eq. (12.101) above. At NLO we find, 
in agreement with Refs. ||9l[T0"1h 
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- 2/9( 165 - 52/x + 276x - 407jc 2 ) - 16(3jc - x 2 )H 2 - 4(24jc - 25jc 2 )H ) 

+ 16(2jc-x 2 )^ 2 -16(1+3x)H ,o + 8(1-2x + 2jc 2 )(H 1 , + H u )]} 

+c 



- l+4x-2(l -2jc + 4x 2 )H -4(1 -2x + 2x 2 )Hi 

4(1 +4x)H + 2/3(3+4/x + 24x-31x 2 ) -4(1 -2x + 2x 2 )Hi } 

8(1+ 4jc)Ho,o + 4/3(3 - 4/x - 12x + 79x 2 )Hi - 4/3 (3 + 96x - 3 1jc 2 )H 
+ (1-2jc + 2x 2 )(8Hi )0 +16Hi,i) -2/3(43 -4/jc + 242jc-281x 2 ) + 16(3jc-x 2 )^ 2 
- 1 6(3jc - x 2 )H 2 ] - 4/3 ( 1 - 8jc + 8jc 2 ) - 4/3 ( 1 - 2x + 2x 2 ) (Ho + Hi) 

8( 1 + 4x)H + 4/3(3 + 4/jc + 2Ax - 3 Ijc 2 ) - 8( 1 - 2x + 2x 2 )Hi 

+ 4/3(1 -2x + 2x 2 )}} + (x,n f + !)/(«/ + 1) . 



+ L Q iC 



dE5) 



(2) 

Here c\ g (x, + 1) is the two-loop gluon coefficient function with massless quarks at /j = Q given 
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+120k* 2 - 144£ 2 2 + 288£ 2 2 *)Ho + 1/9(2379 + 48/x- 4674x + 243(k 2 - 160k 
+320k* - 320k* 2 - 20k - 80k* + 80k* 2 ^ + (1 - 2x)(8H 5 - 16/3H 2 , ( 



in Eq. (B.6) of Ref. IfTTll . Expanding in powers of a s {rir), instead, leads to the additional term 

«S (*, Q\ m 2 ) as{n f +l l4 as{n f ] flg (X; m 2 } _ 2 ^ H (l) (jC) 2 2 ? m 2 } dBl6) 

(3) 

The NNLO contribution //^ g i s new - It i s exact as far as all /^-dependence as well as all terms 
proportional to Lq and are concerned, and is given by 

H^(x,Q 2 ,m 2 ) = 

C F [8/3(1 -2(k + 4x 2 )H , ,o,o + 8/3(1 - 12x + 10x 2 )Hi A0 - 16/3(1 -2x + 4x 2 )H 3j0 

-16/3(1 + 2x-3x 2 )H Ujl - 16/3(1 + 6x -3x 2 )H 2 ,i + 28(2-5x)H 3 
-8/3(3 + 4x)H 2 , +4(5 - 12x)H 4 +4/3(8 - 126x + 20x 2 - 15k + 30k*)H , , 
+4/3(8 - 34x + 24x 2 + 13k -26k* + 26k* 2 )H u 
-8/3(10 - 27x + 16x 2 - 2^ + 4k* - 4^ 2 x 2 )U lfi 

+2/3(44 - 150x - 136x 2 - 16k + 32^ 3 x - 59k + 1 12k* + 32^ 2 x 2 )H ,o 
+4/3(120 - 59x - 12x 2 + 13k - 26£ 2 * + 26k* 2 )H 2 

+1/45(315 - 23670x + 8070x 2 - 1480k + 2960k* - 1440k* 2 - 5315^2 + 5620k* 

3 

0,0) 

+ 1/270(14085 -3 120/x- 129150x+ 104910x 2 -2160k + 4320k* 
+320^ 3 /x - 2 1000k + 45 120k* - 6800k* 2 - 1280k/* - 61215^2 + 38580k* 
+9920k* 2 - 6552^ 2 + 13968^ 2 x- 15696^ 2 x 2 ) + (1 -2x + 2x 2 )(-32/3H 2 , 2 

+ 16/3H 2 ,i ,1 + 16/3H li0 , 0) o - 32/3Hi A0 - 32/3H M)2 + 16/3H UjU ; 

+C F n f [32/3(5 - x - 4x 2 )H 2j0 ,o + 16/3(9 - 6x - 8x 2 )H 3 , + 16/27(19 - 5x - 4x 2 )H Ujl 

-8/9(47 - 4x - 44x 2 )H 0i0 ,o,o - (1 - 2x)(32H , ,o,o,o - 32H 4 , - 32H 3)0 ,o) 
+16/27(19 + 31x-4x 2 )H 2 ,i -4/81(229 -695x + 533x 2 + 108k -216k* + 216k* 2 )H M 
-4/8 1 (283 - 452x + 533x 2 + 108^2 - 216k* + 2 16k* 2 )H 2 
+4/243(355 + 277x - 700x 2 + 432k - 864k* + 864k* 2 ^ 
+4/9(325 + 130x+ 176x 2 )H 2 , -2/27(1697 + 338x+ 1992x 2 )H 0)0 ,o 
+2/81 (6413 + 624 /x - 3688x - 3506x 2 )Hi )0 + 16/27(226 - 24/x - 329x + 1 19x 2 )Hi j0 ,o 
-2/81 (6995 + 6077x- 11582x 2 - 1728k + 3456k*)H , 
-2/243(46904+ 10985x + 46944x 2 - 8208k + 8640k* + 5184k* 2 )H 
-1/14580(7822155 + 612160/x - 5047230x - 2194900x 2 + 933120k - 1866240k* 
-138240k A - 2674080k - 3646080k* - 1244160k* 2 - 378720k + 675360k* 
-733680k* 2 + 12960k 2 + 533952k 2 * - 160704k 2 * 2 ) 

+(1 -2x + 2x 2 )(16/3H 2 , 2 -40/9H 2 , u -80/9Hi,o AO + 16/3H U , 2 -40/9H M)li 
+C A 32/3(1 +4x)H 2 , ,o + 8/3(1 +4x- 5x 2 )H u ,i -4/9(2 + 211x + 72x 2 )H 2 

+8/9(3 - 6x + 23x 2 )H 0i0 ,o - 8/9(6 + 51x + 88x 2 )H 2i0 + 8/9(9 + 8/x + 48x - 65x 2 )Hj A0 
+4/9(19 - 56x)H 3 + 2/8 1 ( 106 - 628/x - 6878x + 68 19x 2 + 1 80k - 360k* + 360k* 2 
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- 1 80^ + 900^* - 900^ 2 )H i +4/27(146-80/*- 1147* +11 28x 2 )H i , 
-4/27 (229 + 592x + 822x 2 + 8l£ 2 + 1 8C 2 *)H ,o - 32/3 (x - x 2 )Hi , 2 
+2 /8 1 (2000 - 4390x + 16303x 2 + 1008^ 3 + 2952^ 3 x + 27^ 2 + 3240^ 2 x + 2448^ 2 x 2 )H 
-1/1215(86835 - 66220/x- 183990x + 151675a: 2 - 5040^ 3 /x + 19440^+27000^* 
+225900^ 3 x 2 - 17280^A + 7965^2 -219780^*+ 104760^ 2 x 2 + 1296^ 2 2 + 49896£ 2 2 x) 
-32/3(x-x 2 )H Uj0 + 64/3(x + x 2 )H_i,_ lj0 - 32/3(x + x 2 )H_ 1 , ,o - 8/3(3x + x 2 )H 2j i 
-4/3(24x + 24x 2 + 5^ 2 + 10^ 2 x + 10^ 2 * 2 )H_i, + 16/3(3^ 3 + 6^ 3 x + 6^ 3 x 2 + 2^ 2 x + 2^ 2 x 2 )H_! 
-32/3(^2 + 2t; 2 x + 2^ 2 x 2 )H_ 1 ,_ 1 -(l + 2x)(16/3H ,o,o,o-32/3H 3 ,o) + 8/3(1 -2x)H 4 
-(l+2x + 2x 2 )(64/3H_ 1 ,_ 1 ,_ 1> o-32/3H_ 1 ,_ 1 ,o,o-16/3H_ 1A o,o-32/3H_ lj 2,o) 
+(1 -2x + 2x 2 )(32/3H 1 , 2j o + 16/3Hi A1 - 16/3H UA0 + 32/3H u , 2 + 32/3H Ujl , 
-16/3H M)M ) + 4/27(67 - 143x + 206x 2 - 45^ 2 + 90^ 2 x - 90£ 2 x 2 )Hi, 



+C A n f 



80/27(2-x + x 2 )Hi j2 +4/27(3-52/x- 165x + 217x 2 +6^2- 12^2X+12^ 2 x 2 )H 



t,o 



+8/27(3 + 23x-6x 2 )H 2 + 16/27(5 -7x + 7x 2 )H Uj0 - 16/9(1 - 12x)H 0: 



0,0,0,0 



+8 /9(4x + 9x 2 )H 3 + 4/8 1 ( 130 + 7x - 88x 2 + 90^ 2 - 180£ 2 x + 1 80^ 2 x 2 )Hi 
+4/243(253 - 890x - 5484x 2 + 432^ 3 - 864^ 3 x - 486^ 2 x 2 )H - 64/3xH 3 



-32/3(1 +4x)H 2i o, + 16/9(5 + 4x + 4x 2 )H_ 1 _ 1 ,o - 8/27(5 + 16x+ l6x 2 )H_ lA0 
-8/9(10 + 42x- 9x 2 )H 2 , + 4/27(16 + 89x + 88x 2 )H , ,o- 8/27(29 - 10x + x 2 )H u ,i 
-8/27(32 + 24/x + 146x - 197x 2 )H lt0) o - 8 /9(6^ 3 + 12^ 3 x + 12^ 3 x 2 - 5£ 2 - 4^ 2 x - 4^ 2 x 2 )H_ 1 
+8/81(56 + 235x + 226x 2 + 27^ 2 + 54^ 2 x + 54^ 2 x 2 )H_ 1,0 + 16/9(2x - 9x 2 )H_ 2 ,o 

13,0 

+4/243 (275 + 36/x + 920x-902x 2 -972^ 3 + 1944^ 3 x- 1944^ 3 x 2 - 90^ 2 - 36£ 2 x 
+36^ 2 * 2 )Hi -2/3645(28335 - 8180/x+ 139290x- 115210x 2 - 8640^ 3 /x + 25920^ 3 
+72090^ 3 x+ 115560^ 3 x 2 + 1620£ 2 - 29880^ 2 x- 3240^ 2 x 2 -23328£ 2 2 -45360£ 2 2 x) 
+4/81(86 + 462x+1147x 2 )H ,o + 4/9(x + 24x 2 )H 2 ,i+ 8/3(^2 + 2^2* + 2^2* 2 )H_i,_i 
+(1 +2x + 2x 2 )(16/3H_ 1 ,_ w , + 8/3H-1 -1,0,0 + 16/3HU -2,0 - 56/9H_i, ,o,o - 16/3H_i A0 ) 

1,1,0,0 

64/9H lil , 2 -8/9H liUi o+40/9H Uil)1 ; 1 



+(1 -2x + 2x 2 )(-8/3H 1A o,o + 16/9Hi )3 -40/9Hi, 2)0 -40/9Hi, 2 ,i +56/9H 



+C| 



8(1 - 2x + 8x 2 )H 0j0 ,o,o,o - 16(1 - 2x - 4x 2 )H 2 , OA o - 32(1 - 2x + 3x 2 )H 3 



3,1,0 

-8(1 -2x + 4x 2 )H 4 , -24(1 -2x + 4x 2 )H 3 , u +4(1 + 16x - 16x 2 )H u , 2 
+4(1 +4x-44x 2 + 25^ 2 -50^ 2 x + 36^ 2 x 2 )H 2j0 + 16(3 - 17x+ 13x 2 )Hi AO ,o 
+8(3- 16x+14x 2 )Hi ;2 ,o + 4(3- 13x + 40x 2 )H 4 + 4(3- 12x + 28x 2 )H 3 , 
-8(3 - 6x + 4x 2 )H 3A o - 32(3 - 6x + 5x 2 )H 2j2 , - 16(3 - 6x+ 10x 2 )H 3j2 
+4(3 - 6x+ 16x 2 )H 5 + 8(3 + 4x - 8x 2 )H 3J +4(3 + 8x - 16x 2 )H 1A1 
+4(3 + 8x - 16x 2 )H ljU , - 4(3 + 16x - 52x 2 )H , ,o,o - 2(5 + 2x + 36x 2 + 3£ 2 
-14£ 2 x + 8£ 2 x 2 )H OA o - 4(5 + 8x - 40x 2 )H 2;0 ,o - 4(7 - 30x + 16x 2 )H 2 , 2 
-8(7- 14x+12x 2 )H 2jl ,i, -8(7- 14x + 20x 2 )H 2jl , u +4(7 + 18x- 18x 2 -4^ 2 
+8^>x - 8^ 2 x 2 )Hi ;0 ,o + 2(7 + 86x - 88x 2 + 36^ 2 - 72^ 2 x + 72^ 2 x 2 )H u , 
+2(7 + 86x - 88x 2 + 60^ 2 - 120^ 2 x + 120^ 2 x 2 )Hi, 2 - 4(9 - 30x + 16x 2 )H 2i i, 
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-8(9 - 18x + 8x 2 )H 2 , 3 + 8(11 - 4x - 31x 2 + 18^ 2 - 36^ 2 x + 40^ 2 x 2 )H 2jl 
-4(1 l-52x + 72x 2 )H 2 , u -8(13- 26x + 28x 2 )H 2 ,i, 2 + 4(15- 56x + 40x 2 )H li3 
+4(15 -56x + 40x 2 )H M , ,o- 8(15 -30x + 28x 2 )H 2 , 2jl -4(19-76x + 72x 2 )H M;M 
+4(31 + 28x - 62x 2 + 40^ 2 - 80^ 2 x + 80^ 2 x 2 )H ljl , 1 + 16(^2 + 2^ 2 x + 2^x 2 )H_ lj _ 1 , 
+2/3(75 -468x + 576x 2 - 88^ 3 + 176^- 176^ - 81^ 2 + 252^ 2 x- 192^ 2 x 2 )H 1 , 
+1 /3(84 - 915x + 816x 2 + 16^ 3 - 32^ 3 x - 32^ 3 x 2 + 120£ 2 - 492^ 2 x + 384^ 2 x 2 )H 2 
- (37 + lOx + 176x 2 - 68^ 2 + 136^ 2 x - 240^ 2 x 2 )H 3 - 1/5(130 + 910x - 1360x 2 + 60^ 3 
+400^ 3 x - 640^ 3 x 2 - 145^2 - 190^ 2 x - 610^ 2 x 2 + 228^ 2 2 - 376^ 2 2 x + 1136£ 2 V)H 
+1/3(189 - 588x+ 1152x 2 - 80^ 3 + 160^ 3 x- 352^ 3 x 2 - 24^ 2 + 222^ 2 x- 312^ 2 x 2 )H , 
+2/3(222 - 387x + 408x 2 - 8^ 3 + 16^ 3 x - 16^ 3 x 2 - 216^ 2 x + 192^ 2 x 2 )H u 
+1/15(585- 4695x + 4080x 2 - 320^ 3 - 1 120^ 3 x + 1920^ 3 x 2 - 300^ 2 - 2460^ 2 x + 1 830^ 2 x 2 
- 1584^ 2 + 3 168^ 2 2 x - 3 168£ 2 V)Hj - 1 /60(23 10 + 2520x + 2340x 2 + 240^ 5 - 480£ 5 x 
+13440^ 5 x 2 -520^ 3 + 1360^*- 10080^ 3 * 2 -5625^ 2 -7650£>x + 3480^ 2 x 2 + 11520^ 2 ln2 
-23040^ 2 ln2x + 23040^ 2 ln2x 2 - 160(^3 + 2240^ 3 * + 3200^^ + 48^ 2 2 - 3792^ 2 2 x 
+8448^ V) - 8(^2 + 2^2* + 2^ 2 x 2 )H_ 1A o - 8(^2 + 2^ 2 x + 4^x 2 )H_2,o 
-8(2^ + 5^ + 3^ 2 )H_ lj0 + 8(C 2 2 + 2^ 2 2 x + 2^ 2 2 x 2 )H_ 1 - 40(1 - 2x)H 2jl A o 
+(1 - 2x + 2x 2 )(32H 1AO ,o,o + 32H li4 - 16H 1A0 - 161^,2,0,0 - 80Hi, 2)2 
— 48Hi,2,i,o — 48Hi ; 2,i,i + 32Hi ; i ; o,o,o — 32Hi ^3 — SOH+i^o — 1 1211^1^,1 
— 112Hijj ;2 — 48Hi j i i i ) i > o — 80Hijjjj' 

+C A C F [8(1 -4x-8x 2 )H_ 2 _i j0 - 8/3(1 - 1 13x + 48x 2 )H 2 , u -4(1-8*- 12x 2 )H_ 2A0 

+96(1 + x + x 2 )H_ 1 ,_ 1 ,_ 1 , - 32(1 +x + x 2 )H 5 + 64(1 + x + x 2 )H 2 , 2 , 

-8(1 +2x + 4x 2 )H_ 2 ,_i A o - 16(1 +2x + 4x 2 )H_ 2 ,_ 1 , 2 - 16(1 +2x + 4x 2 )H_ 2 ,_ 2 , 

+16(1 + 2x + 4x 2 )H_ 2 , OA o + 8(1 +2x + 4x 2 )H_ 2 ,3 + 16(1 + 2x + 4x 2 )H_ 2 , 2 , 

-16(1 +2x + 4x 2 )H_ 3j _ 1 ,o + 8(1 +2x + 4x 2 )H_ 3 ,o,o + 96(1 + 3x)H 2AO ,o 

+64(1 + 4x - x 2 )H 2 , 1A o + 24(1 + 6bc)H 2) i,i,o + 8(1 + 8x + 4x 2 )H_ 1 , 2 , 

+8 /3 (2 + 8/x + 26x - 43x 2 )Hi ;2j o + 8(2 + x - x 2 - 6^ 2 - 12^ 2 x - 12^ 2 x 2 )H_ u 

+32(2 + 3x + 2x 2 )H 3A0 + 4/3(3 + 16/x + 48x-94x 2 )Hi, Uj0 - 8(3 - 18x + 4x 2 )H 2jl , u 

+4/3(3 + 32/x + 264x - 302x 2 )H u - 4/3(3 - 42x - 38x 2 )H 4 

+8(3 - 6x + 8x 2 )H 2j i i2 + 8(3 - 2x + 4x 2 )H 3 , 2 + 8(3 + 6x + 4x 2 )H 3 , u 

+4(4 + 14 X - 2x 2 - 5^ 2 - 10^ 2 x - 20£ 2 x 2 )H_ 2 , - 8(7 + 12x + 8x 2 )H_!,_i, 2 

-8(5 - x - 3x 2 - 3^ 2 - 6^ 2 x- 6^ 2 x 2 )H_ 1 ,_ lj0 + 16(5 + 2x + 6x 2 )H 2 , 2a 

-8(7 + 12x + 8x 2 )H_i,_ 2i0 + 8(7 + 14x + 4x 2 )H 3 , 1 ,o + 4(7 + 16x+12x 2 )H_ 1 ,3 

+4(7 - 4x 2 - 15^2 - 30^ 2 x - 30^ 2 x 2 )H_! A0 - 8/3(8 - 31x + 43x 2 )H 3 , 

-4/3(9 - 156x - 82x 2 )H 3J + 4/3( 13 - 92x + 74x 2 )H 1; i j2 + 8( 13 + 22x + 4x 2 )H 2j3 

+2/3(13 + 178x+16x 2 )H ,o A o-4(19 + 24x + 20x 2 )H_ 1 _i AO + 4(19 + 52x- 106x 2 )H 2 + 

-4(21 + 40x + 35x 2 + 2^ + 4^ 3 x + 4^ 3 x 2 + 13^ 2 + 27^ 2 x + 17^ 2 x 2 )H_i, 

_4(24 + 4/x - 94x + 57x 2 + 48^ 2 - 96^ 2 x + 96£ 2 x 2 )Hi, u + 32(x + x 2 )H_i, 2 j 

+4/3(25 + 40/x + 334x - 404x 2 )H 1A o,o + 8/3(56 + 4/x - 3 lx - 12x 2 )H M , u 
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-2/3(47 + 722* + 480x 2 + 168^ + 96^ + 144^ 2 )H 2;0 + 32(x + x 2 )H- 2 , 2 
+2/9(51 + 1164x - 1868x 2 - 180^ 2 + 360^2* - 432^ 2 x 2 )H 3 + 32(x + x 2 )H_ 3j0 
-2/9(62 - 845x - 1 159x 2 + 246£ 3 + 84£ 3 x + 168£ 3 x 2 + 648^ 2 - 288^ - 492^ 2 )H 2 
+4 /3 (7 1 + 242* - 430x 2 )H 2)0 ,o + 4/3 (91 + 1 12x - 150x 2 )H 2j2 
-4/9(168 - 1317jc+ 1691x 2 + 108^ 2 -432^ + 360^ 2 )H 2jl 
+4/9(297 - 8/jc— 1152x + 818x 2 - 324^ 2 + 648^ 2 x- 648^ 2 x 2 )H 1)2 
+4/9(303 — 8/jc— 1068x + 728x 2 - 180^2 + 360^ 2 x- 360^ 2 x 2 )H M , 

- 1 /9(346 + 2987x - 10622x 2 - 480^ 3 - 624^ 3 x - 288^ 3 x 2 - 147^ 2 - 318^ 2 x + 384^ 2 x 2 )H 0)0 

-1/9(1131 + 1752x + 7624x 2 - 288^2 - 360^ 2 x - 288^ 2 x 2 )H 0j o,o 

+2/9(1269 - 40/x - 1392x+ 124x 2 + 18^ 2 - 36^ 2 x + 36^ 2 x 2 )H! A o 

-2/27(1401 -740/x+ 12504x- 12769x 2 - 1044^ 3 +2088^ 3 x-2088^ 3 x 2 +432^ 2 /x 

+423^ 2 + 2178^ 2 x - 2826£ 2 x 2 )H 1 ,o - 2/27(2241 + 142/x - 3096x + 2369x 2 - 252^ 3 

+504^ 3 x- 504^ 3 x 2 - 144£ 2 /x+ 1026^ 2 - 5508^ 2 x + 5220^ 2 x 2 )H ljl 

+ 1/1620(22990 - 83880/x + 986660x - 545880x 2 - 3240^ 5 + 563760^* - 686880^* 2 

+17820^ 3 - 2689200^ 3 x + 734400^ 3 x 2 - 64800^ 2 /x + 57735^ 2 - 496350^ 2 x - 495000^ 2 x 2 

+ 155520^ 2 In2-311040^ 2 ln2x + 311040^ 2 In 2x 2 - 135000^ 2 ^3 -689040^ 2 ^ 

+254880^2^ 2 + 21816^ 2 2 - 714960^2 2 ^ + 508464^ 2 V) + 8/3(8/x + 30x-47x 2 )H 1 , 2 ,i 

-2/405(56370- 1825/x-220125x+ 151945x 2 +2880£,/x-2925^ 3 +46620^ 3 x 

-3240^2 A - 48960£ 3 x 2 + 26685^ - 164070^ + 123840^ 2 - 35154^ 2 2 + 70308^ 2 2 * 

-70308^2 V)Hi - 1/405(87495 - 171555x + 47570x 2 + 11790^ 3 + 12060^ 3 x + 164520^ 3 x 2 

+15615^2 + 63405^ - 53100^ 2 + V 128^ 2 2 + 66744^ 2 *)H 

+8/3(16/x+ 138x- 151x 2 )H UA0 -4(^ 3 + 2^ 3 x + 4^ 3 x 2 - fc + 16^x 2 )H_ 2 

-8(7^ + 14£ 3 x+ 14^ 3 x 2 - 13^2 - 14^ 2 )H_!,_i +48(^ 2 +2^+2^)11-1-1,-1 

-2(43^ + 56£ 3 x + 52^ 3 x 2 + 18^2 + - - 18^ 2 2 - 36£ 2 2 x - 36fcV)H_i 

2\tt , -ititr , ->r „ , at „2ntt off , or „ , a? 



+32(^2 + 2^ 2 x + 2^ 2 )H_ lj _2 + 16(^2 + 2&C + 4^)H_ 2 ,-i - 8(^2 + 2£ 2 x + 4^ 2 )H_ 3 
-8H 0; o,o,o,o - 16(1 +2x)H 4; i + 16(2 + 5x + 3x 2 )H_ 1A0 , 

+(1 + 2x + 2x 2 )(96H_ 1i _ 1j _ 1i _ 1 ,o-48H_ 1 ,_ 1i _ 1 ,o,o-32H_ 1j _ 1 ,o,o,o-16H_ 1 ,_ 1 ,3 
— 64H_i i _i ) 2,o — 32H_i i _i j 2,i — 16H_i i _2,o,o ~~ 32H_i i _2,2 — 32H_i i _ 3) o 
+32H_ li0 ,o,o,o + 32H_ 1)4 + 48H_i j3 ,o + + 32H_ li2 ,o,o + 32H_i A2 

+3ZH_i A i )0 ) + (1 - 2x + 2x 2 ) (- 16Hi, ,o,o,o - 32H M + 32H lj3 , + 16H lj2 ,o,o 
+96H 1)2i2 + 64H 1 , 2 j,o + 32Hi, 2)1) i - 48H M AO ,o + 32H U , 3 + 1 12H MA0 

+128H MA i - 16H liM)0)0 + 144H U ,!,2 + 80H M , U ,<,) 
,2 



8(l-6x + 6x 2 )H 2 , M -8/3(2 + 8/x+110x-97x 2 )H 2>0 ,o 

1,0,0,0 



- 8/3 (5 - 8/x - 94x + 108x 2 ^ 2 - 4/3 (5 + 8/x + 46x + 60x 2 )H 



-8/3(5 + 8/x + 70x + 84x 2 )H_i A0 + 16/3(5 + 8/x + 94x+ 108x 2 )H_! -1,-1,0 
-8/3(5 + 8/x+106x+120x 2 )H- 1 - 1 , ,o + 8/3(6 + 60x + 73x 2 )H 2 ,2 
-8/3(8 - 8/x- 82x + 93x 2 )H u , 2 - 8/3(8 - 8/x- 82x + 93x 2 )H, iMi0 
-8/3(9 + 8/x + 36x-53x 2 )H u - 8/3(9 + 8/x + 48x-65x 2 )H 1A0 , 
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-8/3(11-20/*- 229x + 24 lx 2 - 21^ 2 + 42^ 2 x - 42^ 2 x 2 )H lj2 + 8/3(6 + 60x + 73x 2 )H 2i i i0 
-8/3 (1 1 - 20/* - 229x + 24 lx 2 + 3^ 2 - 6^ 2 x + 6^ 2 x 2 )Hi , i ,o - 8/3 (3 + 12x + 23x 2 )H 4 
-4/3(11 - 8/x-94x+ 102x 2 )Hi ;2 ,i - 4/3(13 + 24/x + 142x- 190x 2 )H u , , 
+2/9(15 - 4/x- 1440x+ 1321x 2 + 144^ 2 - 288^ 2 x + 288^ 2 x 2 )H Mi i 

-4/3(19 + 8/x+ 166x- 204x 2 )H UjU - 8/3(23 - 38x)H 3 , +4/3(23 - 26x+ 18x 2 )H 0;0 ,o,o 
-4/9(24 - 40/* - 654x - 7 1 8x 2 - 8 1 ^ 2 - 1 62^ 2 x - 162^ 2 x 2 )H_ i >0) o 
+8 /9(24 - 40/x - 546x - 610x 2 + 45^ 2 + 90^ 2 x + 90^ 2 x 2 )H_ 1} _ 1)0 

-2/3(35 + 12x - 225x 2 + 60^ 2 - 96^ 2 x)H 0;0 ,o + 2/27(39 - 134/x - 3 108x + 3536x 2 - 936£ 3 

+1872^ 3 x- 1872^ 3 x 2 - 144^ 2 /x+ 198^ 2 - 1692^ 2 x+ 1836^x 2 )H u 

+8 /9(42 + 39x + 274x 2 + 63^ 2 + 144£ 2 x)H 3 - 4/27(450 + 1008x + 6446x 2 

-360^ 3 - 1440^ 3 x - 36^ 2 /x+ 108^ 2 + 1458^ 2 x + 2349^ 2 x 2 )H 2 - 64(x + x 2 )H_ 1j _ 1j2 

-64(x + x 2 )H_!,_ 2i0 + 32(x + x 2 )H_i, 3 - 2/9(84 + 564x - 533x 2 + 72^ 2 + 288^ 2 x)H 2jl 

+4/9(285 + 80/x + 2052x + 1560x 2 + 90^ 2 + 360£ 2 x)H 2)0 - 4/27(300 - 188/x - 4398x 

-4886x 2 + 216^ 3 +432^ 3 x + 432^ 3 x 2 -108^ 2 /x- 18^ 2 -414^ 2 x-603^ 2 x 2 )H_i,o 

+2 /9(477 - 308 /x - 1788x + 1619x 2 - 18^ 2 + 36^ 2 x - 36^ 2 x 2 )H lj0 ,o 

+2/27(738 - 6573x + 5248x 2 + 720^ 3 - 4464^ 3 x + 423^ 2 - 54^ 2 x + 342^ 2 x 2 )H , 

+2 /27( 1286 - 1465 fx - 18101*+ 18046x 2 + 480^ 3 /x + 246^ 3 + 41 88^ 3 x - 5244^ 3 x 2 

-726^ 2 /x + 348^ 2 -9156^ 2 x + 9333^ 2 x 2 - 1242^ 2 2 + 2484^ 2 2 x-2484^2 2 x 2 )Hi 

-4/27(1575- 1388/x- 16335x+16040x 2 + 180^ 3 -360^ 3 x + 360^ 3 x 2 - 180^2 A - 189^ 2 

-864£ 2 x+ 1233^ 2 )Hi, - 1/405 (163770 + 26240/x+ 1002420x 

+2259620x 2 + 1440^ 3 /x + 48420^ 3 - 61200^ 3 x- 37800^ 3 x 2 + 7200^ 2 A 

+52425^2 + 83880^ + 291330^* 2 + 12312^ 2 2 - 80352^ 2 2 x)H - 32(2x + 3x 2 )H_ 2 ,-i,o 

+ 1/1215(830790 - 666070/x - 11788380x+ 11588020x 2 - 165240^ 5 + 1321920^ 5 x 

+40320^ 3 /x + 31 8330^ 3 + 2261520£ 3 x + 14727 60^ 3 x 2 - 97920^ A + 206550^2 

-459405^ + 2589660^ 2 + 134460^^3 + 450360^2^ + 49248^ 2 2 A-58320^2 2 

+1084104^ 2 x + 55728^2 V) + 48(2x + 2x 2 + ^ 2 + 2£ 2 x + 2^ 2 )H_ ii2 

+ 16(2x + 3x 2 )H_ 2 ,o,o - 8(2x + 3x 2 )H 3J - 32(£ 2 - 4£ 2 x)H_ 2 _i 

-4 /45(360^ 3 /x + 225^ 3 + 4410^ 3 x + 5040^ 3 x 2 + 200£ 2 fx - 120^ 2 + 38 10^ 2 * + 4130^ 2 
+423^ 2 2 + 846^> 2 * + 846^2 V)H_i + 16(3£ 3 - VX&t - 2C, 2 x - 3^ 2 )H_ 2 
+8/3(60^ 3 + 120£ 3 x + 120^ 2 + 8^A + 5^ 2 + H8^+ 132^ 2 )H_!,_i 
+8(24x + 4x 2 - 3^ + 12^ 2 x)H_2,o - 160(^ 2 + 2£* + 2^ 2 )H_ 1j _ 1j _ 1 - 320xH 3 , ,o 
+32(^2 + 2C.2X + 2C 2 * 2 )H_i,_ 2 + ( 1 - 4jc)(-64H_ 2 ,-i,-i,o + 32H_ 2 ,-i ,o,o + 16H_ 2 ,o,o,o 
+32H_ 2 ,2,o - 16Ho,o,o,o,o + 32H 4 , ) + (1 + 4x)(-32H 2 ,o,o,o - 32H 2j3 + 32H 2 , 2 ,o 
+16H 2)2i1 -48H2,i A o+32H 2 , 1)2 + 32H2 > i,i j o-16H 2) i,i,i) 
+(1 + 2x)(16H 5 - 32H 3 , 2 - 32H 3 ,i, ) + (1 +2x + 2x 2 )(-192H_ 1 ,_ 1 ,_ lj _ 1 , 
+128H_i j _i ) _i i o,o + 64H_i j _i ) _i i2 + 64H_i_i_2 i o — 16H-i,-i,o,o,o ~~ 32H_i i _i i3 
+32H_i j _i i 2,o + 64H_i i _ 2 _i ; o — 32H_i j _2,o,o — 16H-i,o,o,o,o ~~ 16H_i 5 4 — 32H_i i3j o 
-32H_ 1 , 2 ,o,o - 32H_ lj2 , 2 - 32H_ 1A1)0 ) + (1 - 2x + 2x 2 )(-32H 1;3 ,o - 16H 1)3) i 
-16Hi, 2A o - 80H lj2)2 - 8QHi A i )0 - 32H li2 , u + 16H 1)1AO ,o - 16H Mi3 - 96H MA0 
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-64H Mi2 ,i 
+ Lq {C F 



8OH11 1,2 — 80Hijjj ; o + 8OH1111 1 

,.2n 



+ 1/9(8^3- I6C3JC+ 16£ 3 x 2 



32/3(5 - 6x + 4x 2 )H_ 2 ,o - 32(5 - I5x- 2x 2 )H ,o,o - 8 /3(33 - 60a: - 4x 2 )H 3 

-8/9(131 - 8/jc- 253x + 168x 2 )H M - 8/9(140- 8/jc- 283a; + 198^)Hi, 
+ 1 6/45 (225 + 1 /x 2 - 20 /x + 40x - 1 55a: 2 + 36a 3 )H_ 1 , - 4/9 (609 - 3 1 2x + 320a: 2 )H 2 
- 8 /45 (2975 - 2660a: + 560a: 2 + 72x 3 - 270^ 2 + 540^ 2 x)H , - 32/3 (£ 2 + 2^ 2 x + ^ 2 x 2 )H_ 1 
-2/27(5247 + 256/x - 8448a: + 3080a: 2 - 216£> +432^>x - 576^ 2 a: 2 )H 1 
-4/45(9213 + 4/jc- 5632x + 1444x 2 - 360^ 3 + 720^3*- 990^ 2 + 2520^ 2 x+ 120^ 2 x 2 )H 
- 1 /405(567861 - 1 1744/x - 828996x + 269954x 2 - 47520^ + 90720£ 3 x - 17280^ 3 x 2 
+2880^ 2 /x- 109620^ 2 +50400^ 2 x- 57600^ 2 x 2 - 5184£>x 3 - 648^ 2 2 + 1296^ 2 2 x) 
-32/3(2 - 4x + 5x 2 )Hi j0 ,o + ( 1 - 2x)(-96H 0j0 ,o,o - 48H 4 - 16H 3;0 - 16H 3J ) 
-8/3(13 -8x + 4x 2 )(H 2i o + H 2j1 ) + (1+2x + x 2 )(-64/3H_ 1 ,_ 1 ,o + 32/3R 



1,0,0 J 



+(l-2x + 2x 2 )(-80/3Hi, 2 - 16H 



1.1,0 



16H 



1,1,1. 



+C F n f 



32/3(5 - 6x + 4x 2 )H_ 2 ,o - 8/3(33 - 60x - 4x 2 )H 3 - 4(35 - 108x - 16x 2 )H , , ( 



-8/9(131 -8/x-253x+ 168x 2 )H M -8/9(140- 8/x-283x+ 198x 2 )H lj0 

+ 1 6/45 (225 + 1 /x 2 - 20/x + 40x - 155x 2 + 36x 3 )H_ 1 ,0 - 4/9 (609 - 3 1 2x + 320x 2 )H 2 



-2/27(5247 + 256/x - 8448x + 3080x 2 - 216^ 2 + 432^ 2 x - 576^ 2 * 2 )H 1 
-4/45(5320 - 3745x + 1 120x 2 + 144x 3 - 540^ 2 + 1080^ 2 x)H , 

-4/45(7593 + 4/x - 5407x + 1984x 2 - 360^ 3 + 720^ 3 x - 990^ 2 + 2520^ 2 x + 120^ 2 x 2 )H 
-2 /405 (229863 - 6952 /x - 3 1 85 13x + 91 102x 2 - 23760^ 3 + 45360£ 3 x - 8640^ 3 x 2 
+ 1440^ 2 /x - 548 10^2 + 25200^2* - 28800^ 2 * 2 - 2592^ 2 x 3 - 324^ 2 2 + 648^ 2 2 *) 
-32/3(^2 + 2^ + ^x 2 )H_! - (1 -2x + 2x 2 )(80/3H lj2 + 16H M , + 16H Ujl ) 
-32/3 (2 - 4x + 5x 2 )Hi j0 ,o + ( 1 - 2x) (- 88H 0j0 ,o,o - 48H 4 - 1 6H 3 , - 16H 3 ,i ) 
-(l+2x + x 2 )(64/3H_ 1 ,_i,o-32/3H_i i0 ,o)- 8/3(13 -8x + 4x 2 )(H 2; o + H 2 j 



-Ca 



16/3 (2 + 4x + 5x 2 )H 1 ,o - 8/3 ( 1 - 66x + 35x 2 )H 2 - 16/3 (2 - 4x + 3x 2 )Hi , ,< 



- 16/9 (4 + 2/x - 7x - 17x 2 )H_ 1 , + 8/9( 10 + 296x - 63x 2 )H 0)0 + 32 /3 (3x - x 2 )H 2 ,o 

-8/9(13 - 4/x- 80x + 87x 2 )Hi ;0 - 8/9(13 -4/x - 80x + 87x 2 )H u + 32/3(3x-x 2 )H 2 ,i 

+16/27(31 + 832x - 325x 2 - 90& + 18^ 2 )H + 32/3 (5x - x 2 )H 3 

-8/27(65 - 20/x - 964x+ 1032x 2 - 9^ 2 + 18£ 2 x)Hi + 32/3x 2 H_ 2 ,o 

+2/8 1 (969 + 412/x + 23556x - 25657x 2 + 324^ 3 - 2376^ 3 x + 648^ 3 x 2 - 144^/x 

+ 108^2 -6624^ 2 x + 37 80^ 2 ) - 8/3(^ 2 + 2^ 2 x + 4^ 2 x 2 )H_i +224/3xH 0j0 ,o 

+16/3(1 +2x)H_! _i, + 16/3(1 + 2x + 2x 2 )H_!, 2 - (1 - 2x + 2x 2 )(32/3H u , + 16/3H MJ ) 



+C A n f 8/3(1 + 26x)H ,o,o - 8/3(1 - 65x + 35x 2 )H 2 - 16/3(2 - 4x + 3x 2 )Hi A o 

+ 16/3(2 + 4x + 5x 2 )H_ 1A0 - 16/9(4 + 2/x-7x-17x 2 )H_ 1 , + 32/3(3x-x 2 )H 2) i 
-8/9(13 - 4/x - 80x + 87x 2 )H M +4/9(39 + 542x - 126x 2 )H , + 32/3 (3x - x 2 )H 2 , 
-4/27(83 -40/x- 1807x+ 1934x 2 - 18^ 2 + 36(; 2 jc)Hi - 8/9(13 -4/x- 80x + 87x 2 )Hi, 
+4/27(89 + 2822x - 1520x 2 - 360^ 2 x + 72^ 2 x 2 )H + 32/3(5x - x 2 )H 3 
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+2/9(181 - 16/x + 2056x-2306x 2 +36£ 3 -264^ 3 x + 72^ 3 x 2 - 16^A+ 12^ 2 -724^ 2 x 
+420^* 2 ) + 32/3x 2 H_ 2 ,o - 8/3 + 2^* + 4^ 2 )H_ i + 1 6/3 ( 1 + 2x)H_ i _ j , 
+16/3(1 +2x + 2x 2 )H_i, 2 + (1 -2x + 2x 2 )(-32/3H u , - 16/3H ljU ) 

144(2 - 8x + 7x 2 )H ua - 32(1 + 2x + 4x 2 )H_ 2 ,2 + 80(3 - 6x + 8x 2 )H 2 , 2 

+96(3 - 6x + 8x 2 )H 2;M - 128(3 + 5x + 2x 2 )H_i A0 - 8(3 + 40x - 50x 2 )Hi A0 
-64(3 + 8x 2 )H_ 3 ,o + 64(4 - 13x + llx 2 )H u , + 16(5 - 26x + 20x 2 )H 2A0 
+32(5 - 10x+ 18x 2 )H 3j0 + 32(5 - 6x + 8x 2 )H_ 2 ,-i,o - 32(5 + 2x + 14x 2 )H_ 2 ,o,o 
+8(5 + 6x + 40x 2 )H OAO ,o + 32(7 - 14x+ 18x 2 )H 2 j,o + 16(8 - 39x + 44x 2 )H 2 , 
+8(9 -39x+100x 2 )H 3 - 32(14 + 2x + x 2 )H_ 2 ,o + 32(14+ 15x + x 2 )H_!,_i, 
+8(15 - 22x + 72x 2 )H 4 + 8(33 - 108x + 88x 2 )H!, 2 + 8(19 - 115* + 126x 2 )H 2 ,i 
+8(16 - lOlx + 92x 2 - 20^ 2 + 72^ 2 x - 64^ 2 x 2 )H 2 + 8(25 - 50x + 88x 2 )H 3jl 
+8(29 - 95x + 55x 2 - 32^ 2 + 64^ 2 x - 64^ 2 x 2 )H lj0 - 128(^2 + 2^>x + 2^>x 2 )H_i,_i 
- (85 + 216x - 272x 2 + 240^ 3 - 480^ 3 x + 480^ 3 x 2 + 40^ 2 - 624^ 2 x + 688^ 2 * 2 )Hi 
+2(155 -546x + 368x 2 - 128^ 2 + 256^ 2 x- 256^ 2 x 2 )H M + 8(29x + 82x 2 )H 0A0 
+2/15(308 - 16/x+ 1593x+ 1464x 2 - 1080^ 3 + 3600^ 3 x - 5280^ 3 x 2 - 540^ 2 
+1860^ 2 x- 6000^ 2 x 2 )H + 16(7^ 3 + 14^ 3 x+ 14^ 3 x 2 + 18^ 2 + 25^ 2 x + 7^ 2 x 2 )H_i 
-1/15(383 - 32/x- 1893x+ 1092x 2 +2760^ 3 - 6480^ 3 x+ 15360^ 3 x 2 + 1920^ 2 
+5320£ 2 x+ 11040^ 2 x 2 - 1152^ 2 x 3 -684^ 2 2 + 4152^ 2 2 x-2304^ 2 2 x 2 ) 
+2/15(585 + 5510x + 4620x 2 - 576x 3 - 900^ 2 + 1320^ 2 x - 4320^ 2 x 2 )H , 
-16/15(915 - 2/x 2 + 1090x+ 165x 2 -72x 3 - 30^ 2 - 60^ 2 x-60^ 2 x 2 )H_i, 
+16(7^ 2 - 2£>x+ 16^ 2 )H_ 2 - 32(2 + 5x + 3x 2 )H_ 1,2 + (1 +2x + 2x 2 )(224H_ 1 ,_ 1A0 
-128HU _i,_i, + 64H_ 1) _i )2 + 128ELi _2, - 64H_ 1AO ,o - 32H_ U ) 
+(1 - 2x + 2x 2 )(- 128Hi,_ 2 ,o + 96Hi A o,o + 256H li3 + 288Hi, 2 , + 352H u ,i 
+160H U A0 + 320Hi, 1)2 + 288H 1 , Ui0 + 384Hi,!,i,i ] 

+C A C F [l6(l + 50x + 20x 2 )H_ 2 ,-i,o - 32(3 + 20x - 6x 2 )H 3 , - 8(15 + 62x - 12x 2 )H 4 

-128(3 + 6x + 4x 2 )H_i _j _i ;0 - 64(3 + 6x + 5x 2 )H_ u - 16(1 + 46x - 16x 2 )H 2j2 

-16(3 + 40x- 16x 2 )H 3a - 16(5 + 18x)H OAO ,o + 32(7 -2x + 4x 2 )H_ 3j0 

-32(7 + 2x + 10x 2 )H_ 2 ,2 - 16( 1 1 + 22x + 18x 2 )H_! , , ,o - 8(13 + 50x + 32x 2 )H_ 2 ,o,o 

+8(26- 8/x- 82x + 93x 2 )H Ujl + 16/3(44- 16/x- 199x + 221x 2 )Hi, 2 

+ 128(2 + 4x + 3x 2 )H_i,_!, 2 - 16/15(55 -2/x 2 -40/x-405x-600x 2 -72x 3 )H_i,_!, 

+8/9(65 + 32/x - 346x + 518x 2 +72^ - 144^ 2 x)H M +4/3(65 - 556x+ 1124x 2 )H 2 , 

+8/3(69 - 32/x - 366x + 398x 2 )H M , +4/3(107 - 556x+ 1006x 2 )H 2jl 

+8/3(161 - 16/x -364x + 303x 2 )Hi A0 + 16/15(340 - 175x-600x 2 -72x 3 )H_ 2 ,o 

+4/9(223 + 64/x - 1520x + 1834x 2 + 288^2 - 576^ 2 x + 288£ 2 x 2 )Hi, 

+8/15(225 -610x + 2305x 2 + 144x 3 )H 0A0 + 4/15(255 - 1660x + 6060x 2 +288x 3 )H 3 

-4/15(295 + 8/x 2 + 80/x+ 1000x+ 1140x 2 +288x 3 )H_ 1A0 

_4/45 (602- 24/x + 10977x- 16274x 2 - 2016x 3 - 1350^-4860^+ 1080£ 2 x 2 )H , 
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+C 2 



+2/45 ( 173 1 + 48/x + 1 656x + 29968x 2 + 540^ 2 + 7560^ 2 x - 21 60^ 2 x 2 )H 2 
- 8 / 1 5 (6 1 5 + 4/x 2 + 670x + 300x 2 + 1 44x 3 )H_ i , 2 + 2 /45 (3445 + 64/x + 1 1758x + 242x 2 
-1080^-20880^- 1440^ 2 - 1530^ 2 + 5760^ 2 x-36360^ 2 x 2 -3456^ 2 x 3 )H 
+8/45(4851 -28/x 2 - 12/x + 6071x + 558x 2 - 1008x 3 + 1890£> + 3780^ + 3060^ 2 )H_ 
+ 1/45 (10349 + 544/x + 22856x-32084x 2 + 17640^ -35280^ + 25200^ 2 
+48t; 2 /x 2 + 2880^ 2 /x- 11880^ 2 + 38040^ 2 x- 38640^ 2 x 2 - 1728^ 2 x 3 )Hi 
+1/45(18123 + 2512/x + 39737x - 67432x 2 - 1920£ 3 A - 2940^ 3 - 38760^ + 360^ 
-8640^ 3 * 3 - 192^ 2 /x - 3462^ 2 + 45256^ - 59936^ 2 x 2 - 8064^ 2 x 3 + 1404^ 2 2 
-8424^ 2 2 x + 6768^ 2 2 x 2 )- 192(3x-x 2 )H 2 , ljl -64(12x-5x 2 )H 2 j j0 
+8/15(765^ + 1530^ + 11 10^ 3 x 2 + 6^ 2 /x 2 + 40^ 2 /x + 560^ 2 + 1075^ 2 x + 900^ 2 x 2 
+216^x 3 )H_i - 64(7^ + 14^+ 10^ 2 )H_i,_i + 8(29^2 + 58^x + 60^ 2 )H_ 2 
+16(1-2x4- 6x 2 ) (Hi )0 ,o,o +H 1)3 ) + (3 + 6x + 4x 2 )(96H_ 1 _ 1)0)0 + 64H_ 1) _ 2 , ) 
+ 1 6 (7 - 14x + 1 8x 2 ) (H 2 ,o,o + Hi,i ,o,o ) - 64( 1 + 2x + 2x 2 ) (H_ i , 2 ,o + H_ i , 2 , i ) 

+(1-2*4- 2x 2 )(96H 1 ,_ 2;0 + 160H 1)2j0 + 96Hi A i + 128H Mi2 + 160H u , lj0 ) 

32(1- 10x-4x 2 )H_ 2i2 - 16(1 +40x-12x 2 )H 2j i,o + 64(1 +2x + 3x 2 )H_i,_i j2 

-16(1 + 40x - 12x 2 )H 2j2 - 32(3 + 28x - 2x 2 )H 3j0 - 16/3(3 + 8/x + 366x - 195jc 2 )H 2 ,o 
+32(3 + 6x + 2x 2 )H_ 1 ,_ 1 ,_i, - 16(1 +2x-2x 2 )H_!_ 2i o + 32(4- 8x + 7x 2 )H UA o 
+32(5 - 10x + 9x 2 )Hi, 3 + 16(5 - 6x + 8x 2 )H_ 2 ,_i j0 + 16(5 + 10x+ l^U-i _ lj0)0 
- 16(5 + 64x - 8x 2 )H 3 j - 8/3(6 + 20/x + 828x - 449x 2 )H 2; i - 4(9 + 266x - 24x 2 )H 4 
- 8/3 (7 + 28/x + 3 17x + 389x 2 )H_ i , ,o + 8/3 ( 15 - 28/x - 384x + 433x 2 )Hi, 2 
+8(9 - 62x - 16x 2 )H_ 2j0 ,o - 8(9 + 18x + 22x 2 )H_ li0 ,o,o + 8(13 - 26x + 22x 2 )H lj0 ,o,o 
+8/9( 10 - 58/x - 3023x + 3243x 2 - 216£ 2 + 432^ 2 x - 360^ 2 x 2 )Hi j0 + 16(11-5 lx)H 0;0 ,o,o 
+8/3(13- 16/x- 166x-218x 2 )H_ lj2 + 8(13- 18x)H_ 3j0 - 16(7 + 14x+16x 2 )H_i, 3 
+8/3(19 -8/x + 92x+100x 2 )H_i _ U )- 4/3(21 +2047x-708x 2 )H 3 
+8/3(29 -24/x-382x + 424x 2 )H Mi i + 8/3(37- 28/x- 428x + 477x 2 )H M , 
-4(37 + 186x- 8x 2 )H 2j0 ,o - 8/3(57 - 8/x + 123x + 260x 2 )H_ 2;0 

^ T !0,0,l 

+8/9(128 + 240/x- 257x- 272x 2 +90^ 2 + 180^ 2 x + 252^ 2 x 2 )H_ 1 , 
+8/9(405 + 104/x-3975x + 4600x 2 +63^2 + 774^ 2 x- 144^ 2 x 2 )H 2 
+2 /27(469 + 3996 /x - 54734x + 5269 lx 2 - 3780^ 3 + 7560^ - 6048^ 3 * 2 + 1 152^/x 
-198^ + 12168^- 13788^ 2 )Hi -32(l+2x + 2x 2 )(H_ lj2 ,o + H_i,2,i) 
+2/9(3849- 14584x+ 18688x 2 + 162^2 + 4140^ 2 x-432^ 2 x 2 )H , 
-2/27(7837 + 1 136/x + 68257x - 37209x 2 - 3672^ 3 - 13824^ 3 x - 288^ 2 /x - 378^ 2 
-32418^ 2 x + 12744^ 2 x 2 )H - 192(3x-x 2 )H 2 , u - 16(^ 2 + 2^ 2 x + lO^jH-i-i 
+4/3(18^3+ 36^ 3 x + 120^ 3 x 2 + 24^ 2 /x - 7^ 2 + 424^ 2 x + 536^ 2 x 2 )H_ i 
+ 1 /405 (448905 - 286600 fx - 2646840x + 25 19 1 85x 2 + 5 1 840^ 3 A + 597240^3* 
-8100^ - 67 6080^^ 2 +48960^ 2 A - 145800^ 2 + 1338480£ 2 x- 1656000^ 



1.0 



-4/3(61 +72/x + 63 lx - 822x 2 )Hi i0 , -4/3(83 + 1648x- 198^)^ 
-8 /9(79 - 71 A - 3443x + 366 lx 2 - 180^ 2 + 360^ 2 x - 288^ 2 )Hi 
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* 2 : 

-2/27(87 + 40/x + 894x- 1048a; 2 - 144^ + 72^ 2 ) - 8/3(1 + 8x)H 
- 1 6/3 (3x - x 2 )H 2 - 8/3 ( 1 + 2x + 2x 2 )H ! , + 8/3 ( 1 - 2x + 2x 2 ) (Hi, + Hi, 



+25434^ 2 2 - 2916^ 2 2 x + 57Q24& V) + (1 - 2x+ 2a: 2 ) (96H 1)2)0 - 48Hi _ 2 ,o + 1 12Hi A i 
+ 112H u , 2 + 112Hi, u ,o + 96H u , u ) + 8(C2 + 34^ + 24^ 2 )H_ 2 ]} 
+ L 2 ,{c>(« / +1) 4/3(31 - 56x-8x 2 )H , +4/9(112 -8/x-215x+130x 2 )H! 

+1 /27(4977 + 256/x - 7632x + 2156x 2 - 216^ 3 + 432^ 3 x - 468^ 2 + 288^ 2 x - 144^ 2 x 2 ) 
+2 /9(57 1 - 236x + 244x 2 - 36^ 2 + 72^ 2 x)H + 4/3( 13 - 8x + 4x 2 )H 2 
+(l-2x)(24H ,o,o + 8H 3 ) + (l-2x + 2x 2 )(32/3H 1 ,o + 8H M ) 

+C A (n f + 1) [4/9(13 - 4/jc- 80x + 87x 2 )Hi - 4/9(13 + 136x- 49x 2 )H 

Io,o 

16(1 + 2x + 4x 2 )H_ 2 ,o -64(1-2* + 3x 2 )H 2 ,i + 16(2 + 5x + 3x 2 )H_i, 

- 1 /2(33 + 92x - 104x 2 - 56^ 3 - 16^ 3 x - 256^ 3 x 2 - 48^ 2 + 88^ 2 x - 320^ 2 x 2 ) 

-32(2 - 7x + 5x 2 )H u - 16(3 - 6x + 10x 2 )H 2 , - 4(3 + 2x + 24x 2 )H , ,o 

-8(5 - 18x+ 10x 2 )H!, -4(6 - 31x + 40x 2 )H 2 -4(9 - 18x + 40x 2 )H 3 

-(1 +2x + 2x 2 )(32H_ 1 ,_ 1 , - 16H_i, ,o) -2(11 -6x+ 18x 2 - 18^ 2 + 20^ 2 x- 80^ 2 x 2 )H 

-(65 - 126x + 36x 2 - 64^ 2 + 128^ 2 x- 128^ 2 x 2 )Hi -4(9x + 32x 2 )H , 

-16(^2 + 2^ 2 x + 2^ 2 x 2 )H_i - (1 - 2x + 2x 2 )(32H[ , 0;0 + 80Hi, 2 + 80H U , + 96H MJ ) 
+C A C F [24(1 +4x)H ,o,o - 8(1 -26x+ 12x 2 )H 2 ,i + 8(1 +22x- 8x 2 )H 2 , 

+8(3 + 20x - 8x 2 )H 3 - 8/3(11 - 76x+ 157x 2 )H 2 - 8/3(14 - 49x+ 141x 2 )H , 
-8/3(30-8/x- 123x+131x 2 )H 1 , 1 -4/3(53- 16/x-232x + 248x 2 )Hi, 
+ 1/9(109 + 1 126x - 2320x 2 - 216£ 2 - 1440^ 2 x + 576^ 2 x 2 )H 
-2/9(191 + 32/x - 1078x + 1052x 2 - 216^ 2 + 432^ 2 x - 432^ 2 x 2 )H 1 
-1/6(193 -48/x-536x + 52x 2 - 192^ 3 x 2 - 176^ 2 + 1360^ 2 x-2512^ 2 x 2 ) 
-16(^ 2 + 2^ 2 x + 2^x 2 )H_! + (1 + 2x + 2x 2 )(- 12H_ 1)0 + 8H_i,o,o + 16H_ li2 ) 
- (1 - 2x + 2x 2 ) (40Hi A o + 48Hi, 2 + 48H U , + 48H Ujl ) 

32(1 + 7x-x 2 )H 2 , - 4/3(1 - 16/x- 118x- 146x 2 )H_i, - 32(1 -4x)H_ 2 ,o 

+16(l + 16x-4x 2 )H 2 ,i - 16(3 -16x)H , ,o- 4/3(17 -24/x-310x + 352x 2 )H!,i 
-4/3(5 - 24/x- 238x + 280x 2 )H!, + 8/3(9 + 8/x + 255x- 126x 2 )H 2 + 16(3 + 22x)H 3 
+4/3(23 + 436x - 54x 2 )H , - 2/9(1009 + 104/x - 3512x + 5183x 2 +216^2 + 1008£ 2 x)H 
-2/9(59- 168/x-4204x + 4549x 2 - 144^ 2 + 288^ 2 x-288^ 2 x 2 )Hi 
+1 /9(509 - 1260/x + 7318x - 6639x 2 - 432^ 3 + 144^ 3 x -216^2 - 4704£>x + 3024^ 2 x 2 ) 
-32(^2 + 2^2X + 2^ 2 x 2 )H_i-(l+2x + 2x 2 )(32H_ 1 ,_i,o-24H_i,o,o-16H_ 1 ,2) 
-(l-2x + 2x 2 )(24H 1 ,o,o+48H 1 ,2 + 48H u ,o+48Hi, u )]} 

" L Q { C f ("/ + 1 ) - 4/9 ( 1 1 - 4x)H - 4/27 ( 84 - 8 /x - 147x + 62x 2 ) - 8/3 ( 1 - 2x)H , 
-16/9(1 -2x + 2x 2 )Hi 



+C 2 
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+C A (n f + l) 8/9(1 +4x)H + 4/27(3+4/x + 24x-31x 2 ) -8/9(1 -2x + 2x 2 )Hi 
+C| [8/3(1 -4x)Hi + 16/3(1 -2x + 4x 2 )H 2 + 4/3(1 -2x + 8x 2 )H , 
+1 /3(1 1 -2k- 16^2 + - 64^ 2 x 2 ) - 4xH + (1 - 2x + 2x 2 )(16/3H!,o + 32/3Hi i ) 



+C A C F 



46/45(15 - 10x + 195x 2 +72x 3 )H 0j0 - 16/3(3 - 20x + 21x 2 )H 2 - 32/3(1 - 2x + 4x 2 )H , 
46/3(6 -14x+13x 2 )Hi, - 64/3(1 -2x + 3x 2 )H 3 + 32/45(90+1 /x 2 + 40x + 36x 3 )H_ 



4/9(11 -4(k+ 106x 2 )H - 8/3(1 + 4x)H ,o - 8/3(1 + 10x-4x 2 )H 2 
-2/3(5 + 16x- 10x 2 -4^ -40^+ 16^ 2 ) + (1 -2x + 2x 2 )(16/3H lj0 + 32/3H u ) 
+4/9(19 - 8/x- 104x + 106x 2 )H! 

16/3(1 - 8x)H , -4/9(1 + 16/*+ 118x- 146x 2 )Hi - 32/3(1 + 4x)H 2 
-4/9(11 + 8/*+ 152x- 18x 2 )H + 2/27(369 -272/jc- 1980x+ 1883x 2 + 144^> + 576£ 2 x) 
+(1 - 2x + 2x 2 )(16/3H!, + 32/3Hi,i)l } 

- L M |c F 128/3(1 +x 2 )H_ 2 ,o - 8/3(11 - 44x + 42x 2 )H M - 16/3(5 - 10x+ 12x 2 )H 2 j 

2\TT 

0,0,0 
1,0 

4/3(17 - 72x + 54x 2 - 8^ 2 + 16^ 2 x- 32^ 2 x z )Hi - 64/3(^ 2 + 2^ 2 x + ^ 2 )H_i 
4/45(401 + 8/x- 294x+ 1098x 2 - 240^ + 480^ - 720^ 2 )H - 64/3x 2 H lj0 ,o 
4/45(497 _ 8/x-462x- 132x 2 -420^ 3 - 120^ 3 x- 1200^ 3 x 2 - 180^> + 880^ 2 x 
1260^ 2 * 2 - 288^ 3 ) - (1 +2x + x 2 )(128/3H_ lj _ 1 ,o - 64/3H_i )0 ,o) 
(1 - 2x + 2x 2 )(64/3H 2 , + 32Hi j2 + 64/3H u , + 32H ljU ; 

32/3(5 - 6x + 4x 2 )H_ 2 ,o - 16/3(7 - 5x + 2x 2 )H 2 ,i - 16/3(11 -4x- 4x 2 )H 2 , 

-8/3(29- 148x- 16x 2 )H 0j o,o - 8/3(31 -56x- 12x 2 )H 3 - 8/3(93 -47x + 39x 2 )H 2 
-8/9(180- 16/jc- 312x+163x 2 )H!, -4/9(195 - 16/x-402x + 250x 2 )H M 
+ 1 6/45 (225 + 1 /x 2 - 20/x + 40x - 155x 2 + 36x 3 )H_ i ,o - 32/3 (^ 2 + 2^ 2 x + ^ 2 x 2 )H_ i 
-8/45(1635 - 2215x - 495x 2 + 72x 3 - 270^ 2 + 540^ 2 x)H , - 32/3x 2 Hi j0 ,o 
-4/27(2379+ 128/x-3837x+1312x 2 -36^2 + 72^ 2 x- 144^ 2 x 2 )Hi 
-8/135(7677 + 6/x- 10773x + 5176x 2 - 1395£ 2 + 3600£ 2 x + 540£ 2 x 2 )H 

- ( 1 - 2x) (48H 0) o,o,o + 48H 4 + 32H 3 , + 1 6H 3J ) 

-(1 +2x + x 2 )(64/3H_ 1 ,_ 1 ,o - 32/3H_i i0 , ) - (1 -2x + 2x 2 )(16Hi j2 + 32/3H Uj0 + 16H!, U ; 
-4/135(24233 - 272/x- 36683x+ 12497x 2 - 2250^ 3 +7380^ 3 x- 2880^ 3 x 2 + 240^>/x 

-8370^2 + 3750^ 2 x - 3510^ 2 x 2 - 432^ 2 x 3 + 594^ 2 2 - 1 188£ 2 2 x) 

+C A [64 /3 (1 + 3x)H , ,o + 32/3 ( 1 + 8x - 2x 2 )H 3 + 32/3 ( 1 + lOx - 2x 2 )H 2 , 

+16/3(3 + 6x + 8x 2 )H_i A0 - 16/3(3 - 6x + 4x 2 )H!, ,o + 8/9(3 + 16/jc + 168x - 199x 2 )Hi, 
-8/3(4 -72x + 73x 2 )H 2 - 8/9(9- 8/x -96x+ 107x 2 )H u -32/9(9 + 2/x-3x- 13x 2 )H_i j( 
+4/27(123 - 104/x+ 1380x- 1588x 2 +72^ 2 - 144^ 2 x + 72^ 2 x 2 )H 1 +64/3x 2 H_ 2 ,o 
+4/27(393 + 1584x - 2980x 2 - 72^ 2 - 576^ 2 x + 144^ 2 x 2 )H + 64/3(3x -x 2 )H 2 ,i 
+8/3 (92x- 57x 2 )H ,o - 64/3x 2 H_ lj _ 1 , - 32/3(1 +2x + 3x 2 )^ 2 H_i 



+C F rif 
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+32/3(1 +2x + 2x 2 )H_ 1i2 -(1-2x + 2x 2 )(32/3H 1; 2+64/3H 1j1i o)+4/27(188 + 40/x 
+490x - 799;t 2 + 1 80^ + 72^ 3 x + 216^ 2 - 48£ 2 A + 72 ^2 - 1224^ 2 x + 1314^ 2 jc 2 ) 
+C A n f \ 1 6/3 ( 1 + 8x)H 0i o, + 1 6/3 ( 1 + 10* - 2x 2 )H 2;0 - 1 6/9 (9 + 2 /jc - 3x - 1 3x 2 )H i , 

-4/9(17 - 16 jx - 208x + 239x 2 )H lj0 + 8/3(3 + 6x + 8x 2 )H_i A0 - 4/3(2 - 128* + 69x 2 )H 2 
-8/3(3 - 6x + 4x 2 )Hi A0 + 4/9(26 + 476x- 119x 2 )H , -4/9(49 - 8/x- 176x + 187x 2 )H M 
-4/27(54 - 40/x - 1743x + 1870x 2 - 36^ 2 + 72^ 2 x - 36^ 2 x 2 )H! - 32/3x 2 H_ 1 ,_i, 
+4/27(87 + 2553x - 1462x 2 - 360^ 2 x + 72^ 2 x 2 )H + 32/3(3x - x 2 )H 2J + 32/3x 2 H_ 2;0 
+32/3(5x-x 2 )H 3 - 16/3(^ 2 + 2^ 2 x + 3^ 2 x 2 )H_! + 16/3(1 +2x + 2x 2 )H_ 1)2 
-(1-2* + 2x 2 )(16/3H 1 , 2 + 32/3H u , ) +4/27(212 - 26/x + 2950x - 3163x 2 + 126^3 
-36£ 3 *+ 108^ 3 x 2 -24^ 2 /x+ 18^2 - 1116£ 2 x + 62l£ 2 x 2 ) 



+C A C F 



48(1 - 12x + 8x 2 )H 3;1 + 144(1 -6x + 4x 2 )H 2;U +64(1 -2x + 3x 2 )H u 



,1,0,0 

-32(1 +4x)H 0;0 ,o,o- 64(1 +4x + 2x 2 )H_ 2 ,o,o + 32/15(5 - 20x- 255x 2 -36x 3 )H_ 2 ,o 

+16(5 - 34x + 20x 2 )H 2 , lj0 - 128(1 +x 2 )H_ 2 , 2 - 32(1 + lOx- 4x 2 )H 4 

+16(3 - 42x + 20x 2 )H 2j2 + 16(19 - 4/x - lOlx + 105x 2 )H U;1 + 32(3 - 6x + 8x 2 )H li3 

+ 16/3(8 - 166x + 147x 2 )H 2;0 + 32/3(12 - 2/x - 33x + 36x 2 )Hi, , 

+28/3(13- 158x+ 144x 2 )H 2jl + 16(14- 4/x- 71x + 76x 2 )H u 

+4/45(26 + 20/x - 2025x + 13734x 2 + 900^ 3 - 1440^ 3 x 2 - 660^ 2 + 7800^ 2 x - 9900^ 2 x 2 
- 1728^ 2 x 3 )H + 4/3(50- 40/x -732x + 765x 2 -24^ 2 + 48^ 2 x- 144^ 2 x 2 )H lj0 
+ 16/15(55 - 690x + 825x 2 +72x 3 )H 3 + 8/15(55 - 330x + 840x 2 + 144x 3 )H 0A0 
+2/3(61 - 48/x - 2524x + 2562x 2 - 192^ 2 + 384^>x - 576^ 2 x 2 )H u 
+8/3(79 - 16/x - 356x + 372x 2 )H u , - 32/15(90 + 1/x 2 + 40x + 36x 3 )H_ u 
+32/15(85 + 1/x 2 + 20/x + 270x + 255x 2 + 36x 3 )H_i,_ lj0 - 192(2x-x 2 )H 3j0 
+4/15(221 + 8/x -3614x + 5973x 2 - 180^ 2 + 1560^ 2 x-960^ 2 x 2 )H 2 + 128H_ 3i o 
+2/15(689 - 196/x - 8019x + 71 16x 2 - 1680^ 3 x 2 + 8^/x 2 + 320^/x - 1000^ 2 + 6360^ 2 x 
-7080^ 2 x 2 -288^ 2 x 3 )Hi - 16/15(175 + 2/x 2 +20/x + 310x + 255x 2 +72x 3 )H_ 1A0 
-4/45(987 - 24/x + 5902x - 8 109x 2 - 1584x 3 - 360^ 2 - 3600£ 2 x + 1440^ 2 x 2 )H , 
+8 /45(1206 - 22/x 2 - 12/x + 2756x + 603x 2 - 792x 3 + 720^ 2 + 1440^ 2 x + 720^ 2 x 2 )H_i, 

-2/45(3587 -980/x-8847x + 7065x 2 - 480^ A + 3030^ -10020^ + 29880^ 2 
+4320^ 3 x 3 + 96^/x + 1326^ 2 - 32708^ 2 x + 35838£ 2 x 2 + 3168£>x 3 - 2124^ 2 2 + 1296£ 2 2 x 

-3168^2V) + 128(4x + x 2 )H_2,-i,o + 64x 2 H 1AO ,o + 64(2^ + 4^ + 3^ 2 )H_ 2 
+16/15(210^3 + 420k* + 210^ 3 x 2 + 3^/x 2 + 20^/x + 265^ 2 + 350^ 2 x + 255^ 2 x 2 
+ 108kx 3 )H_i + (l + 2x + x 2 )(192H_ 1 ,_i i o,o-256H_ 1 ,_i,_ 1 ,o + 128H_i i _ 1 ,2 
+128H_!,_ 2 ,o - 64H_ 1A0 , - 64H_ 1;3 ) + (1 - 2x + 2x 2 )(64H 2A0 + 160H Uj0 
+288Hi, 2 ,i + 256H M , 2 + 224H M i n + 384Hii.ii) - 256(^2 + 2 £ 2 x + ^)H_i,_i 



+C 2 



32(1 - 4x)H_ 3i o + 64(1 - 4x - x 2 )H_ 2 ,2 - 96(1 + 7x - x 2 )H 2 



2,0,0 

128(1 + 9x -x 2 )H 3i i - 96(1 + lOx - 2x 2 )H 2 ,i ;0 - 64(1 + 13x - 3x 2 )H 2 , 2 
64( 1 + 18x - x 2 )H 3j0 - 32(1 + 30x)H 4 + 32(2 - 1 lx - 3x 2 )H_ 2 ,o,o 
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+32(2 - 4x + 3x 2 )Hj AO ,o - 32(2 + 4x + 5x 2 )H_ im o + 32(3 - 16x)H OAO ,o 

f 1,0,0 

,0,0 



-8(3 + 8/x + 266x - 180x 2 )H 2 ,i - 4/3(3 + 64/x + 690x - 804x 2 )Hi 
-32(3 + 6x + 7x 2 )H_ 13 + 32(4 - 8x + 7x 2 )H u + 32(5 - lOx + 9x 2 )H u , 



+8(9-8/x-96x+107x 2 )H Ujl - 16/3(9 - 4/x - 102x- 138x 2 )H_! _ lj0 

+ 16/3(11 -24/x-274x + 316x 2 )H M , - 16/3(15 -4/x+108x + 64x 2 )H_ 2 ,o 

+8/3(17 - 40/x - 478x + 548x 2 )H lj2 - 8/3(20 + 24/x + 926x - 401x 2 )H 2 ,o 

-8/3(23 + 844x- 72x 2 )H ,o,o - 8/3(30+ 1070x- 309x 2 )H 3 

+8 /3 (37 - 8/x - 1 30x - 1 98x 2 )H_ ! , 2 + 4/3 (45 - 32/x - 474x - 642x 2 )H i A o 

+8/9(228 + 116/x- 1209x- 1354x 2 + 126^ 2 + 252^ 2 x + 324^ 2 x 2 )H_ li0 

+2/9(563 - 200/x - 1 1740x + 12621x 2 - 576^ 2 + 1 152^ 2 x - 864^ 2 x 2 )H u 

-2/9(597 + 80/x + 13813x - 33051x 2 - 576^ 3 + 144^ 3 x - 96^ 2 /x - 360^ 2 - 10248£ 2 x 

+3708£ 2 x 2 )H + 2/9(631 -808/x- 16856x+ 17721x 2 -504^ 2 + 1008^ 2 x-720^ 2 x 2 )H 1 , 

+2/9(858 + 256/x - 16848x + 19733x 2 + 288£> + 3312^ 2 x - 576^ 2 x 2 )H 2 

-2/9(1052- 1168/x+13775x- 14269x 2 +312^ 2 -72^ 3 x-432^ 3 x 2 -4512^ 2 x + 4920^ 2 x 2 

-432^ 2 /x + 36^ 3 )H! + 8/3(45^ 3 + 90^ 3 x + 132^ 3 x 2 + 12^ 2 /x - 46^ 2 + 232^ 2 x + 336^ 2 x 2 )H_! 

-32(2^ 2 - IIC2X- 4^ 2 x 2 )H_ 2 - 32(3^ 2 + 6^ 2 x + lO^^H-i _i - 192(3x - x 2 )H 2 , M 

+2/9(2546 - 14800x + 14385x 2 + 144^ 2 + 3744^ 2 x)H ,o + 64(3x + 2x 2 )H_ 2j _ li0 

+(1 + 2x + 3x 2 )(96H_ 1 ,_ 1i o,o + 64H_ 1 ,_ 1i2 )-(1+2x + 4x 2 )(64H_ 1 ,_ 1i _ 1 ,o-32H_ 1 ,_ 2> o) 

+(1 - 2x + 2x 2 )(160Hi j2 , + 128H 1)2 ,i + 160H M , 2 + 192H Ujl , ) + 1 /45(14385 - 5000/x 

+34270x - 42935x 2 - 4260^ 3 - 15960^ 3 x - 31920^ 3 x 2 + 2080^ 2 /x - 8580^ 2 - 197330^ 2 x 2 

+120120^ 2 x- 1584^ 2 2 - 31968^ 2 2 x + 4320^ 2 2 x 2 ) - (1 +2x + 2x 2 )(32H_ li2 ,o + 64H_i,2,i; 

16/3(1-6*+ 10x 2 )H + 32/3(l -2x + 4x 2 )H ,o + 16/3(3 - 6x + 8x 2 )H 2 
+8/3(7 - 24x + 20x 2 )H! + (1 - 2x + 2x 2 )(64/3H!, + 64/3H M ) 
+4/3(21 -40x+ 14x 2 - 12^2 + 24^ 2 x- 32^ 2 x 2 ) 

16(2 - x)H 2 + 8/3(29 - 52x - 16x 2 )H , + 8/3(91 - 33x + 28x 2 - 6^ 2 + 12^)H 
+4/9(183- 16/x-342x+ 184x 2 )Hi +4/27(2397 + 128/x- 3693x+ 1132x 2 - 108^ 3 -216£ 2 
+216^ 3 x+ 108^ 2 x) + (1 - 2x)(48H OA o + I6H3) + (1 - 2x + 2x 2 )(32/3Hi, + 32/3H u ; 

+C A \l 6/9(3 - 4/x - 60x + 67x 2 )Hj - 64/3(1 + 3x)H , - 16/3 (24x - 25x 2 )H 
-4/27(285- 104/x + 642x-904x 2 - 288^ 2 x+ 144^ 2 x 2 ) - 64/3(3x -x 2 )H 2 
-32/3(1 +2x + 2x 2 )H_i i0 + 32/3(l-2x + 2x 2 )(Hi ; o + H M ) 

16/3(1 + 8x)H 0j0 + 8/9(13 -4/x - 80x + 87x 2 )Hi - 8/9(13 + 136x-49x 2 )H 
-32/3(3x-x 2 )H 2 - 16/3(1 +2x + 2x 2 )H_!,o + 16/3(1 -2x + 2x 2 )(H 1 ,o + H u ) 
-4/27(87 + 40/x + 894x- 1048x 2 - 144^ 2 x + 72^ 2 x 2 ) 

+C A C F [32(1 +4x)H 0A0 - 32(1 - 14x + 8x 2 )H 2J + 16(1 + 16x- 8x 2 )H 3 
+16(1 +22x- 8x 2 )H 2;0 -4(7 - 138x+ 152x 2 )H 2 - 8/3(11 - 82x+ 168x 2 )H , 



L Q < C F 



+C F n 



+C A n f 



50 



-32/3(13 - 4/x -71x + 72x 2 )H u -8/3(31 - 16/x- 212x + 228x 2 )H!, 
-2/3(45 - 32/x - 1060x + 1032x 2 - 192^ 2 + 384^ 2 x - 384^ 2 x 2 )H! 
+4/3(47 + 248x-492x 2 - 12^ 2 - 192^ 2 x + 96^ 2 x 2 )H 
-(1 -2x + 2x 2 )(64H 1A0 + 128Hi, 2 + 128Hi,i, + 19ZHi,i,i) 
+2/3(71 + 24/x + 35 lx - 402x 2 + 72^ 3 + 192^ 3 x 2 + 42£ 2 - 828 ^ 2 x + 912^ 2 x 2 ) 

+C 2 64(1 + 7x -x 2 )H 2j0 - 8/3(1 - 16/x- 118x- 146x 2 )H_i, - 64(1 - 4x)H_ 2j0 

+32(1 + 16x-4x 2 )H 2 j - 32(3 - 16x)H 0;0 ,o - 8/3(5 - 24/x- 238x + 280x 2 )Hi, 
+32(3 + 22x)H 3 + 8 /3(23 + 436x - 54x 2 )H , + 16/3(9 + 8/x + 255x - 126x 2 )H 2 
-8/3(17 - 24/x - 310x + 352x 2 )H M + (l+2x + 2x 2 )(48H_ 1 ,o,o-64H_ 1 ,_i i0 + 32H_ 1 , 2 ) 
-4/9(59 - 168/x - 4204x + 4549x 2 - 144^ 2 + 288^ 2 x - 288^ 2 x 2 )H 1 
-4/9(1009+ 104/x-3512x + 5183x 2 +216^ 2 + 1008£ 2 x)H 

-64(^2 + 2^ 2 x + 2^ 2 x 2 )H_! - (1 - 2x + 2x 2 ) (48Hj , 0)0 + 96H lj2 + 96H Uj0 + 96H Ujl ) 
+2/9(509- 1260/x + 731 8x- 6639x 2 -432^ + 144^ 3 x-216^ 2 -4704^ 2 x + 3024^ 2 x 2 ) } 
+ L 2 , jc F [- 4/3(1 -4x)- 8/3(1 -2x + 4x 2 )H - 16/3(1 -2x + 2x 2 )Hi 
+C F n f [ - 8/3(5 - x - 2x 2 )H - 2/9(165 - 16/x - 282x + 124x 2 ) - 8(1 - 2x)H, 
-8/3(1 -2x + 2x 2 )Hi 

16/3(1 +4x)H + 8/9(3 + 4/x + 24x-31x 2 )- 16/3(1 -2x + 2x 2 )Hi 
8/3(1 +4x)H + 4/9(3 + 4/x + 24x-31x 2 )- 8/3(1 -2x + 2x 2 )Hi 



0.0 



+C A n f 

+c A c F 



8/3(4 -4/x-41x + 42x 2 )Hi -8(1 +4x)H , - 8(1 + 10x-4x 2 )H 2 



+C 2 



+2/3(1 1 - 58x+ 168x 2 )H - 1/3(41 + 52c - 60x 2 - 24^ 2 - 240^ 2 x + 96^ 2 x 2 ) 
+(l-2x + 2x 2 )(16H 1 ,o + 32H u ) 

16(1 - 8x)H 0;0 -32(1 +4x)H 2 -4/3(1 + 16/x+ 118x- 146x 2 )Hi 
-4/3 ( 1 1 + 8/x + 152x - 1 8x 2 )H + 2/9(369 - 272/x - 1980* + 1 883x 2 + 144^ 2 + 576^ 2 x) 
+(l-2x + 2x 2 )(16H 1 ,o + 32H u )]}} 

8/9(1 - 8x + 8x 2 ) - 8/9(1 -2x + 2x 2 )(H + Hi) 

+C F n f U/3(9 - 36x- 8x 2 )H , + 4/3(9 - 8x 2 )H 2 + 4/9(81 - 8/x - 135x + 62x 2 )Hi 
+4/9(108 - 153x - 34x 2 - 18^ 2 + 36£>x)H + (1 - 2x)(8H OA o + 8H 3 ) 
+4/9(197 + 4/x - 482x + 28 lx 2 - 18^ 3 + 36^ 3 x - 27^ 2 + 24^ 2 x 2 ) 



+Ca 



16/9(3 -4/x-72x + 79x 2 )H! - 32/3(1 +4x)H 0j0 - 16/9(3 + 96x- 3 lx 2 )H Q 



- 8 /9 (43 - 4/x + 242x - 28 lx 2 - 72^ 2 x + 24^ 2 x 2 ) - 64/3 (3x - x 2 )H 2 
+(l-2x + 2x 2 )(32/3H 1 , + 64/3H 1 , 1 ; 



+C A n f 4/9(3 - 4/x - 72x + 79x 2 )Hi - 8/3(1 +4x)H , - 4/9(3 + 96x -3lx*)Ho 
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+C 2 A 



-2/9(43 - 4/x + 242x - 281a: 2 - 72^ 2 x + 24^ 2 x 2 ) - 16/3(3* - x 2 )H 2 
+(1 -2x + 2x 2 )(8/3Hi j0 + 16/3H U ) 

16(1 + 16jc)H 3 - 16(1 -8x)H 0j o,o + 8/3(6 + 4/x + 279x- 82x 2 )H 2 

+4/3(11+ 440x - 1 8x 2 )Ho,o - 2/9 (477 + 24 fx - 4932x + 25 1 3x 2 + 72^ 2 + 1 1 52^ 2 x)H 
+2/9(8 1 + 224/x + 5412x - 5783x 2 + 216^ 2 - 432^ 2 x + 432^ 2 x 2 )n 1 
+(11 - 16/jc- 214x + 242x 2 )(-4/3H 1 , - 8/3H u ) + (1 + lOx - 2x 2 )(16H 2 , + 32H 2;1 ) 
- 1 /9(259 + 284/x - 17962x + 17419x 2 + 144^ 3 + 2304^ 3 * + 96^ 2 /x + 144^ 2 + 6696£ 2 x 
- 1968^ 2 x 2 ) - ( 1 - 2x + 2x 2 ) ( 16Hi j0 ,o + 48Hi, 2 + 48Hi j, + 96H U j ) 

4/3(9 - 8x 2 )H - 4/9(81 - 8/x- 135x + 62x 2 ) - 8(1 - 2x)H , 



f L Q \ C F n f 



+C A 



32/3(1 +4x)H + 16/9(3 + 4/x + 24x-31x 2 ) -32/3(1 -2x + 2x 2 )H! 



+C A n f 



8/3(1 +4x)H + 4/9(3 + 4/x + 24x-31x 2 )- 8/3(1 -2x + 2x 2 )Hj 



+Cj 



16(1 -8x)H , -4/3(1 + 16/x+118x- 146x 2 )Hi -4/3(11 + 8/x+ 152x- 18x 2 )H 
+(1 - 2x + 2x 2 )(16Hi,o + 32Hi,i) +2/9(369 - 272/x - 1980x + 1883x 2 + 144^ 2 + 576^ 2 x) 
-32(l+4x)H 2 ] +8/9(1 -2x + 2x 2 )}} +fl^°(jc) + c^ g (x,n f + l)/(n f + l). ©7) 

Again, c^(x,n^ + 1) denotes the gluon coefficient function with massless quarks at fj = Q, now 
at three loops, as determined in Eq. (B.9) of Ref. [17], where all contributions proportional to flf x 
are omitted as discussed in Section 13.31 The ^-independent part of the massive OME has been 

approximated in Eqs. (13.491 ) and (13.501 ) and is denoted by Decoupling of the heavy quark in 



a s in the MS scheme leads to the additional matching terms H60H61I 

flg (x, Q\m 2 ) -\h, H® (x, Q\ m 2 ) - { (f C A - % C 
+ (fc A + 2C F ) L,-\ L 2 } flg (x, 2 2 , m 2 ) . 



F ^ 

dj8) 



(£) ps (1) ps 

The corresponding expansion coefficients H 2 ^ at LO and NLO are given by = and 



<< ps (x,2 2 ,m 2 ) 



C F [4(1 +x) (4^3 - Ho,o,o + 2U 2fi ) + 4/3(3 + 4/x - 3x - 4x 2 )H 1;0 
+ 2/3(3 + 1 5x + 8x 2 ) H , - 2/27 (9 + 1 1 2/x + 279x - 400x 2 ) - 4/9 (2 1 + 33x + 56x 2 )H 
+ L Q C F [8( 1 +x) (£ 2 - 2H ,o - H 2 ) - 4/3(3 + 4/x - 3x - 4x 2 )Hi 
- 8/9(39- 13/x-30x + 4x 2 ) + 16x 2 H 

4( 1 +x)H + 2/3 (3 + 4/x - 3x - 4x 2 ) 
8(1 +x) (C2 - H ,o - H 2 ) - 4/3(3 + 4/x - 3x - 4x 2 )H! 



+ LqC f 



+ L u iC 
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Sx 2 



-4/3(3 + 15x-4x 2 )H - 8/3(10 -l/x-x 
+ L Q C F 8(1+x)H + 4/3(3+4/x-3x-4x 2 ) } + 



,(2),ps 

■2,q 



x,n f +i)/(n f +l) 



dj9) 



where c^' ps (x,n f + 1) denotes the two-loop pure-singlet quark coefficient function at /u = Q as 
written down in Eq. (B.7) of Ref. ifTTI 

Finally the new third-order (NNLO) pure-singlet coefficient function for heavy-quark produc- 
tion at Q 2 3> m 2 reads 

H^(x,Q 2 ,m 2 ) = 



+C F n 



C F 8/9(3 + 15* + 8x 2 )H ,o, + 16/9(5 -Ax- 6x 2 )H 2 , + 8 /9(7 + 22/x - Ix - 22x 2 )H M 

+16/9(8 + llx + 2x 2 )H 3 + 8/9(45 + 51x+ 16x 2 )H 2J -8/27(145 + 79x + 66x 2 )H 2 

+8 /27(63 - 20/x- 1 17x + 74x 2 )Hi )0 - 16/27(64 + lOOx + 75x 2 + 9£ 2 + 9^ 2 x)H ,o 

+8 /8 1 (463 - 473x + 784x 2 + 342^ 3 + 342^ 3 x + 27^ 2 + 135^ 2 x + 72^ 2 x 2 )H 

+8/81(93 -254/x + 177x- 16x 2 )Hi -8/243(2400- 1480/jc- 1842x + 922x 2 

- 144 /x^ 3 +999^ +2511^+ 1332^ 3 x 2 - 360/x£ 2 - 1386^ 2 - 1116^ 2 x + 252^ 2 x 2 

-243^ 2 2 - 243^ 2 2 x) + (3 + 4/x - 3x - 4x 2 ) ( 16/9H 1A o - 8/9Hi, 2 - 8/9H Uj0 - 32/9H U j ) 

+16/3(1 +x) (2H 3 , - H OA o,o - H 3) i + 2H 2;0 ,o - H 2 , 2 - H 2) i,o - 4H 2 , U ] 

16/27(5 -Ax- 6x 2 )H 2il +8/27(49 + 13x - 24x 2 )H , ,o - 32/27(5 - 4x - 6x 2 )H 3 

- 16/81 (43 + 60/x + 47x - 150x 2 + 36/x^ 2 + 27^ 2 - 27^ 2 x - 36^ 2 x 2 )Hi 

-16/9(13 + 7x-4x 2 )H 2j0 - 16/3(5 -7x + 2x 2 )H!, + 8/81(63 -20/x - 117x + 74x 2 )Hi > i 

+8/243(769 + 967x-300x 2 - 108^ 3 - 108^+ 180^ 2 - 144^ 2 x- 216^ 2 x 2 )H 

+ 16/729(2775 - 752/x - 3720x + 1697x 2 +432/x^ 3 - 1647^ 3 - 1566^ + 54^ 3 x 2 + 900^ 2 

+576^ 2 x + 999^ 2 x 2 + 729^ 2 2 + 729^ 2 2 x) + 16/8 1 ( 137 + 29x + 54x 2 + 36^ 2 + 36£ 2 x)H ,o 

+(3 + 4/x-3x-4x 2 )(16/9Hi, 2 -16/9Hi i0 ,o + 16/9Hi il ,o-8/27H M ,i) 

- 1 6/8 1 ( 1 00 + 64x + 1 1 lx 2 + 54^ 2 + 54^ 2 x)H 2 + ( 1 + x) ( 80 /9H 0;0 ,o,o - 64 /9H 4 

-32/3H 3 ,o + 32/9H 3 ,i - 32/3H 2j0 ,o + 32/3H 2j2 + 32/3H 2 , lj0 - 16/9H 2 , U ) 
+C F 2 8/3(3 + 9x + 4x 2 )H OA o,o-8/3(3-33x-8x 2 )H 3 , 1 - 16/3(3 -4x 2 )H 4 

+16(4 + 5x)H 2j2 + 8/3(5 - 6/x + 7x - 6x 2 )H Ujl - 8/3(9 - 33x- 8x 2 )H 3i0 

+4/9(6 - 87x - 260x 2 + 18^> + 18^ 2 x)H 2J + 32/3(6 + 9x - 2x 2 )H 2A0 

-8/3(9 + 3x + 4x 2 )H 2 , u + 16/3(9 + 6x-4x 2 )H 2jl , + 4(15 + 8/x- 15x-8x 2 )H!, 2 

+8/3(43 + 10/x-31x- 22x 2 )Hi j0 ,o- 8/9(51 +228x+ 148x 2 + 18£ 2 + 18£ 2 x)H 2 , 

+2/9(93 + 465x - 304x 2 + 144^ 2 + 144^ 2 x)H 3 + 4/3(41 + 24/x - 41x - 24x 2 )H u , 

-4/27(288 + 71/x + 225x - 584x 2 - 36/x^ 2 - 27^ 2 + 27^ 2 x + 36^ 2 x 2 )H M 

-4 /27 (53 1 + 1530x + 604x 2 - 36^ 3 - 36^ 3 * + 35 1 ^ 2 + 405 ^ 2 x + 36^ 2 x 2 )H 2 

-4/27 ( 1 179 - 370/x + 135x - 944x 2 + 72/x^ 2 + 54^ 2 - 54^ 2 x - 72^ 2 x 2 )Hi )0 

-2/27 ( 1548 + 1899x - 88x 2 + 72^ 3 + 72^ 3 x - 324^ 2 - 108^ 2 x + 432^ 2 x 2 )H , 

-2 /27 (28 19 + 753/x - 728x - 2844x 2 - 48 /x^ 3 - 36^ 3 + 36^ 3 * + 48^ 3 * 2 + 396 /x^ 2 + 477 L, 2 
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-801£ 2 * - 72^ 2 x 2 )H 1 - 4/9(207 + 225* + 212x 2 - 45^ 2 - 45^)H 0;0 ,o 
- 1 /27 (4594 - 4030x - 4080x 2 + 648^ 3 - 2304^ 3 x + 96^x 2 + 1251^2 + 4203^* - 1 104^ 2 x 2 
+ 1 080^ 2 + 1 080£ 2 2 *) H o - 2/ 1 35 ( 14 1 80 + 3495 fx - 1 1 1 1 5x - 6560x 2 - 8 1 00^ 5 - 8 1 00£ 5 * 
+7350^ + 32070^* + 5400^ 3 * 2 + 2700/x^ 2 - 6525^ 2 - 15165£ 2 x - 7660k* 2 + 4140^ 3 
+4140^ 3 x + 2160^ 2 2 - 2484^ 2 2 x - 2664^ 2 V) + (3 + 4/x -3x- 4x 2 )(8Hi, 0A0 
+16/3H 1)3 + 16/3Hi A0 + 16/3H u ,i + 16/3H UA o + 16/3H Uil , + 8/3H liMil ) 
+8(1 +x)(2H 3 ,i,i + 6H 2 ,o, ,o -H 5 - 2H 4 , - 2H 4 ,i - 2H 3 , 2 + 2H 3>1 , +4H 2 , 3 + 4H 2 , 2 ,o 
+4H 2j2i i +4H 2jl ,o, +4H 2 , U , + 2H 2 ,i ;U ; 



+C A C F 



32/9(6 - 10/jc - 3x- 19x 2 )H_i _!, - 16/9(6 - 10/x-3x- 19x 2 )H_i A0 



+32(1 - 2x)H 3i o, + 16/9(6 + 10/jc + 3x - 19x 2 )Hi >2 + 16/9(6 + 10/jc + 3x - 19x 2 )H Uj0 

-8/9(15 + l/x+12x-28x 2 )H Uil +8/9(132 -113/x- 141x+122x 2 )Hi A0 

+8/3(17 - 4x)H , OA o - 4/3(36 + 47x + 36x 2 + 18^ 2 - 18^ 2 x)H OA o 

-16/3(17 -jc)H 3 , -4/9(42 + 93x + 20x 2 + 18£ 2 + 18£ 2 x)H 2 ,i 

-8/3(19 + 8/x + 19x)H 2 ,o,o + 8/9(150 + 40/x + 240x + 176x 2 +9^ 2 + 9^ 2 x)H 2; o 

-4/27(228 + 69x+ 1 12x 2 - 180^3 - 180£ 3 x - 36/x^ 2 + 99^ 2 + 234£>x + 36^ 2 x 2 )H 2 

-4/27(246 + 67 /x + 15x - 328x 2 + 36/x^ 2 + 27^ 2 - 27^ 2 x - 36^ 2 x 2 )H u 

-4/27(300- 188/x-240x-728x 2 - 108/x^. + 81^ 2 + 81^ 2 x - 108^ 2 x 2 )H_ !i0 

+4/27(855 -68 lx + 2132x 2 + 468^ 3 - 504^ 3 x + 306^ 2 - 180^ 2 x)H , 

-2 /8 1 (984 + 1771 /x - 255x - 2500x 2 - 720/x^ - 540^ 3 + 540^ 3 x + 720^ 3 x 2 + 414 /x^ 2 

+ 1053^ 2 - 81^2*- 1386^ 2 )Hi + 16/3(4/x^ 2 - 3^ 2 - 3^+4^ 2 x 2 )H_i _i 

-4/27( 1596 - 1348/x - 3324x + 307 6x 2 - 36/x^ 2 - 27^ + 27^ 2 x + 36^ 2 )H li0 

-2/405(85410+ 13120/x + 46695x + 271020x 2 +720/x^ 3 +38700^3 -2340^ + 3600/x^ 

+ 1 3275^2 + 5445^2* + 2 1330^ 2 x 2 + 1 620^ 2 2 - 123 12^ 2 2 x)H 

+2/405(147095 - 108825/x-415850x + 377580x 2 -24300^ 5 + 56700£ 5 x + 6840/x^ 3 
+61830^ + 104220^* + 69120£ 3 x 2 - 15450/x^ + 36900^ - 51570^ 2 x + 55230£>x 2 
+9990^3 + 11610^ 3 x + 8208/x^ 2 2 -7398^ 2 2 + 18360£ 2 2 x+ 1080^ 2 2 x 2 ) 
-8/9(36/x^ 3 -27^ 3 -27^ 3 x + 36^x 2 + 20/x^- 12^ + 6^x + 38^ 2 x 2 )H_i -8xH 3j i 
+24(2x - £ 2 + t,2x)K-2,o ~ 32(^2 - Ca*)H_ 2l _i + 8(^ 2 + k*)H 3 + 48(^ 3 - k*)H_ 2 
+8/3(3 + 4/x - 3x - 4x 2 )(2H lj2 , + H lj2 ,i - H MA0 + 2H M , 2 + 2H ljU;0 - Hujj) 
+16(1 +x) (2H 2j2 ,o + H 2)2 ,i - H 2 ,i A o + 2H 2 ,i, 2 + 2H 2 ,i,i, -H 2 ,i,i,i) 

+8/3(3 -4/x + 3x-4x 2 )(H_ 1 ,o A o-4H_ 1 ,_ 1 ,_ 1 ,o + 2H_ 1 ,_ I ,o,o + 2H_ 1 , 2 ,o) + 8(1 -x)(H 2 , u 
+4H_ 2,-1,0,0 — 8H_ 2 ,-i,- i,o + 2H_ 2,0,0,0 +4H_ 2, 2 ,o — 2Ho A o,o,o +4H4,o) 
+ L e | C F 1 6/9 (29 + 59x - 1 2x 2 )H ,o - 64/9 (3 + 1 /x + 3x + x 2 )H_ i , + 32 ( 1 + x)H A o 
+ 16/27(233 + 155x - 60x 2 )H + 16/81(687 + 124/x - 507x - 304x 2 - 36/x^ 2 - 108^ 2 x) 

+C f .« / [l6/9(8 + llx + 2x 2 )H 2 - 64/9(3 + l/x + 3x + x 2 )H_!, - 8/9(3 + 4/x-3x-4x 2 )H M 

+ 16/81(615 + 68/x- 489x - 194x 2 + 27^ 3 + 27 ^ 3 x - 36/x^ 2 - 72^ 2 - 207 £ 2 x - 18^ 2 x 2 ) 
+ 1 6/27 ( 1 90 + 88x - 79x 2 )H + ( 1 + x) (32H , ,o - 16/3H 2) i ) 
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+c A c F 



+8/27(27 + 20 fx -9x- 38x 2 )Hi + 16/9(29 + 59x - 12x 2 )H 0j0 
+C F 2 64/3(3 + l/x + 3x + x 2 )H_ u -64/3(9- l/x + 9x-x 2 )H_i _ 1)0 

+ 1 6/3 ( 1 5 - 22* + 52x 2 )H , ,o + 32/3(15 + 5* - 2x 2 )H_ 2 , - 1 6/3 (2 1 - 6* - 44x 2 )H 2 ,o 
-16/3(30 + 51x-56x 2 )H 3 - 8/3(57 + 15x-76x 2 )H 2il -64(l/x-x 2 )Hi, 0>0 
-4/3(325 - 20/x- 313x + 8x 2 )H u - 4/3(331 - 32/x- 331x + 32x 2 )H!, - 4/135(606 
+48/x- 32859x + 6298x 2 - 1080^ 3 + 3240^ 3 x - 5400^ 2 - 10980^ 2 x+ 10080^ 2 x 2 )H 
-8/9(723 -333x-220x 2 - 108^ 2 - 108^ 2 x)H 2 + 16/45(1065 + 4/x 2 + 1105x- 36x 3 )H_!, 
-4/27(4455 + 292/x-5733x + 986x 2 -648/x^ 2 + 108^ 2 - 108^ 2 x + 648^ 2 x 2 )Hi 
-2/135(5159 - 4056/x - 20639x + 19536x 2 +2880/x^ 3 -9180^ -39420^ 
-43380^ 2 -6540^ 2 x+ 13200^ 2 x 2 + 864^ 2 x 3 +4968^ 2 2 + 7128^ 2 2 x) + 128H_ 3 ,o 
-2/45(7695 - 1525x - 3680x 2 - 288x 3 - 6480^ 2 - 6480^ 2 x)H , - 64^ 2 (1 - x)H_ 2 
-32/3(l/x+ 15 + 15x + x 2 )^ 2 H_! -8/3(3 + 4/x-3x-4x 2 )(10Hi ;2 + 10H M ,o + 7H M; i) 
+(1 +x)(128H_ 1A o - 168H ,o,o,o - 288H 4 - 240H 3j0 - 224H 3jl - 64H 2 ,o,o 
- 160H 2j2 - 160H 2jli0 - 1 12H 2 , M ) - 64(1 - x) (2H_ 2 ,-i, - H_ 2 , , ) 

24 ( 1 - 4/x + x - 4x 2 )H_ i i0 ,o - 8/3 ( 1 + 20/x + 4x - 28x 2 )H 2 , ! - 32 (2 + 5x)H 3i0 
- 16/3(3 + 8/x + 6x - 14x 2 )H 2j0 + 8(3 - 25x)H 4 - 16(5 + 1 lx)H 3;1 
+ 16(7 - 3x)H_ 3 ,o + 16(11- 12x)H 0j0 ,o,o + 16/9(12 - 31/x - 57x + 76x 2 )Hi j0 
-4 /3 (2 1 + 87x - 40x 2 )H 3 - 8 (34 + 41x)H A o - 8/3 (39 - 8/x - 3x + 40x 2 )H_ 2j0 
-4/3 (69 + 56/x - 69x - 56x 2 )H 1A o + 4/9 (8 1 - 106/x - 261x + 286x 2 )H U 
+4 /9 (90 1 + 208 jx + 73x + 846x 2 +216^2 + 21 6^ 2 x)H 2 

-4/27(909 - 2210/x + 333x + 968x 2 - 360/x£> - 432£ 2 + 432^ 2 x + 360^ 2 x 2 )H 1 
+4/9(1297 - 1040x+ 1508x 2 - 54^ 2 + 342£ 2 x)H , + 8 (8/x + 3 + 3x + 8x 2 )£ 2 H_! 
-4/27(8014 + 568/x + 3922x + 932x 2 - 1620^ 3 - 1512^ 3 x - 144/x^ 2 - 189^2 - 837^ 2 x 
+360^ 2 x 2 )H + 16/9(156+ 118/x+ lllx + 73x 2 )H_i ;0 -48(l -x)^ 2 H_ 2 
+4/135(8980 - 23020/x - 2635x + 16675x 2 +4320/x^ 3 - 1260^ 3 - 3150^ 3 x-6840£ 3 x 2 
+3960/x£ 2 - 13515^ 2 + 5565^ 2 x - 12690^ 2 x 2 + 945^ 2 2 + 1107^ 2 2 x) -64(l/x + x 2 )H_i, 2 
-40/3(3 + 4/x - 3x - 4x 2 ) (H li2 + H u>0 + H Mjl ) + (1 +x) (48H_ 1 ,_ 1 ,o - 136H 2A o 
-80H 2 , 2 - 80H 2 , lj0 - 80H 2jM ) + (1 -jc)(32H_ 2) _i,o + 112H_ 2 , , + 64H_ 2j2 )] } 

+ Lq lc P (n f + 1) - 16/9(8 + llx- 4x 2 )H - 8/9(35 - 2/x - 23x- 10x 2 ) - 32/3(1 +x)H 0j0 

+C F 2 \ - 8/3(6 + 3x + 28x 2 )H 0j0 +4/9(198 - 315x - 4x 2 - 180^ 2 - 180£ 2 x)H 

+2/9(253 + 36/x - 409x+ 120x 2 - 144^ 3 - 144^ 3 x -216^2 - 72^ 2 x + 336^>x 2 ) 
+4/9(369 - 32/x - 333x - 4x 2 )Hi + (3 + 4/x - 3x - 4x 2 ) (8Hi, + 8H U ) 

+ 1 6/3 (9 + 3x - 14x 2 )H 2 + ( 1 + x) (48H 0;0 ,o + 80H 3 + 48H 2;0 + 48H 2 
+C A C F [l6(l + 4x)H 3 + 8/3(3 + 8/x - 8x 2 )H 2 - 8/9(5 1 - 44/x - 69x + 62x / )Hi 
-16(3 -4x)H 00 o - 16/9(82 + 13/x- 53x+ 101x 2 +9^ 2 + 18^ 2 x)H 
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+4/27(1314 - 901 fx - 594x + 181x 2 - 216^ 3 + 108^ 3 x - 54^ 2 - 108^ 2 x + 144£ 2 x 2 ) 
-16/3(3 -4/x + 3x-4x 2 )H_!, + 16/3(3 + 4/x-3x-4x 2 )(H 1 , + H 1 ,i) 

+(1 +x)(224/3H ,o + 32H 2 ,o + 32H 2i i) - 32(1 -jc)H_ 2 ,o] } 
+ L 3 \c F {n f + l) [8/27(3 + 4/x-3x- 4x 2 ) + 16/9(1 +x)H 



+C F 



16/9(3x + 4x 2 )H - 4/9(23 - 23x - 24^ 2 - 24^ 2 x) - 16/9(3 + 4/x - 3x - 4x 2 )Hi 



-16/3(1 +x)(H ,o + 2H 2 ; 



+C A C F 



16/3(1- 2x)H 0;0 - 8/9( 1 1 + 4/x + 8x)H - 8/9(3 + 4/x - 3x - 4x 2 )Hi 



-4/9(45 -44/x-45x + 44x 2 + 12^ 2 + 12^ 2 x)- 16/3(1 +x)H 2 } 

{ C F 32/27 ( 19 + 5/x - 76x + 52x 2 + % 3 + - 6 /x^ 2 - 1 8 ^ 2 x + 1 8 ^ 2 x 2 ) 

-64/9(3 + 1 /x + 3x + x 2 )H_i j0 + 16/27(105 - 99x - 88x 2 - 36^> - 36^ 2 x)H 
+32/27 (39 - 1 3 /x - 30x + 4x 2 )H i + 32/3 ( 1 + x) (2H ,o,o + 2H 3 + 2H 2j0 + H 2 , i ) 

+(3 + 4/x - 3x-4x 2 )(32/9H lj0 + 16/9H U ) + 32/3(5x-2x 2 )H , -64/3x 2 H 2 



-C F n f 



64/9(4 +13x-2x 2 )H , - 64/9(3 + l/x + 3x + x 2 )H_ 1)0 + 16/3(5 -7x + 2x 2 )Hi 



+ 16/9( 13 + 7x - 4x 2 )H 2 + 32/9(22 + 15x - 13x 2 - 3^ 2 - 3^ 2 x)H 
+ 16/9(55 - 39x - 16x 2 + 12^ + 12^ 3 * - 4/x^ 2 - 13£ 2 - 19^ 2 x + 4^ 2 x 2 ) 
+ 1 6/9(3 + 4/x - 3x - 4x 2 ) (Hi ;0 - ! ) + 32/3 ( 1 + x) (2H , ,o + H 3 + H 2;0 - H 2 , 

+C F 2 64/3 (3 + 1 /x + 3x + x 2 )H_ i j2 - 8 /9(345 + 44/x - 255x - 1 34x 2 )Hi, 

- 64/3 (9 - 1 /x + 9x - x 2 )H_ i i ,o - 1 6 (9 + 6x - 8x 2 )H 2j0 + 32/3 ( 1 5 + 5x - 2x 2 )H_ 2 ,o 
- 16/3(15 + 69x - 32x 2 )H 3 - 16/3(36 + 33x - 28x 2 )H 2J + 16/3(3 - 34x + 20x 2 )H OA o 
-8/9(477 + 16/x- 351x- 142x 2 )H M -32/3(l/x+ 15 + 15x + x 2 )^ 2 H_i 
-4/9(799 + 58/x- 637x- 220x 2 - 216/x^ 2 + 36^ 2 - 36^ 2 x + 216^ 2 x 2 )Hi 
-4/9(825 - 333x- 668x 2 - 216^ 2 - 216^ 2 x)H 2 + 16/45(1065 + 4/x 2 + 1105x- 36x 3 )H_!, 
-4/45(1172+ 16/x-5863x- 1244x 2 - 1620^ 3 - 180^ 3 x-900^ 2 -4740^ 2 x+ 1920£ 2 x 2 )H 
-4/45(2089 - 526/x - 3879x + 2316x 2 - 240/x^ - 900^3 - 5760^ + 1080^ 3 * 2 - 4125^ 2 
-2755^ 2 x + 3340^ 2 x 2 + 144^ 2 x 3 + 180^ 2 2 + 540^ 2 2 x) - 64(1 -x)^ 2 H_ 2 + 128H_ 3 ,o 
-4/45(2385 + 1915x-3260x 2 - 144x 3 - 1980^ 2 - 1980^ 2 x)H , - 64(l/x-x 2 )H 1 , , 
- 16/3(3 + 4/x - 3x - 4x 2 ) (5Hi, 2 + 4H U , + 3H u ,i ) - 64(1 - x)(2H_ 2 _i,o - H_ 2 , 0)0 ) 
+(1 +x)(128H_ 1 ,o, - 80Ho,o,o,o - 176H 4 - 160H 3 , - 19ZH 3) i - 64H 2i0 , 
- 160H 2j2 - 128H 2jli0 - 96H 2jM ) 

+C A C F [32(1 - 6x)H 4 -96(1 +2x)H 3a - 8/3(31 +24/x + 46x-24x 2 )H 2i0 

+ 16/3(3 + 8/x + 3x + 8x 2 )H_ 1 ,_i ;0 + 32(3 - 4x)H , ,o,o + 8/3 (7 - 24/x + 1 6x + 32x 2 )H 2 , 1 
+8(3 - 8/x + 3x - 8x 2 )H_!,o, - 16/3(15 - 4/x + 18x + 8x 2 )H_ 2i0 - 32(1 +7x)H 3j0 
-8/3(21 + 16/x- 21x- 16x 2 )Hi A o -32/3(10+ 19x-4x 2 )H 3 - 16/3(28 + 79x)H 



10,0,0 



-8 /9(69 - 50/x - 132x + 1 13x 2 )H U + 8/9(214 + 64/x - 17x + 489x 2 + 90^ 2 + 90^ 2 x)H 2 
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+8/9(117-202/*- 270* + 355x 2 )H lj0 + 8 /9(460 - 347* + 451x 2 - 36^ 2 + 1 80^)H 0i0 
-8/27(1014- 826/x-987x + 799x 2 - 144/x^ 2 - 189£> + 189£ 2 x+ 144£ 2 x 2 )Hi 
- 8/27 ( 1677 + 60/x - 132x - 3263x 2 - 324^ 3 + 108^ 3 x - 72/x£ 2 - 360^ 2 - 360£ 2 x 
+ 144^ 2 x 2 )H + 32/9(57 + 28 /x + 3x - 26x 2 )H_!, + 8(8/x^ 2 - 3^ 2 - 3^ 2 x + 8^ 2 x 2 )H_! 
-4/135(8485 + 3560/x - 53080x + 41035x 2 + 4230£ 3 + 10170^ 3 x + 2160^ 3 x 2 - 1440/x^ 2 
+6420^ 2 - 870^ 2 x+ 14670^ 2 x 2 + 1404^ 2 2 +4536^ 2 2 x) - 80(1 -x)^ 2 H_ 2 
+32/3(3 -4/x + 3x-4x 2 )H_!, 2 - (3 + 4/x-3x-4x 2 )(16H 1 , 2 + 64/3H M ,o + 16H Ui i) 
-16(l+x)(5H 2 ,o ) o + 6H 2)2 + 8H 2) i,o + 6H 2il , 1 ) 

+ 16(1 -jc)(5H_ 2 ,o,o -2H_2,-i,o + 4H_ 2 ,2 + 2H_ 3 ,o) 

h L Q |c f 64/3x 2 H - 32/27(39 - 13/x - 30x + 4x 2 - 9^ 2 - 9^ 2 x) - 16/9(3 + 4/x - 3x - 4x 2 )Hi 
-32/3(1 +x)(2H ,o + H 2 ) 



+C F n f 



+C 



32/9(8 + llx- 4x 2 )H - 16/9(35 - 2/x- 23x - 10x 2 ) - 64/3(1 +x)H, 



o.o 



16(3 + 4/x-3x-4x 2 )(Hi j0 + H 1 , 1 ) +4/9(315 -531x-232x 2 - 288^2 -288^)H 

+4/9(299 + 36/x - 395x + 60x 2 - 72^ 3 - 72^ 3 x - 252^ 2 - 252^ 2 x + 240^ 2 x 2 ) 

-16/3(3 - 6x + 20x 2 )H , + 8/9(333 + 8/x - 225x - 116x 2 )Hi + 16/3(21 + 21x - 20x 2 )H 2 

+16(1 +x) (5H ,o,o + 8H 3 + 6H 2 , + 6H 2 ,i ) 

+C A C F [32(1 +4x)H 3 + 16/3(3 + 8/x - 8x 2 )H 2 - 16/9(51 -44/x- 69x + 62x 2 )Hi 

-32(3 - 4x)H OA o - 32/9(82 + 13/x - 53x + 101x 2 + % 2 + 18^ 2 x)H 

+8 /27( 1314 - 901 fx - 594x + 181x 2 - 216^ 3 + 108^ 3 x - 54^ 2 - 108^ 2 x + 144^ 2 x 2 ) 

-32/3(3 - 4/x + 3x - 4x 2 )H_!, - 64(1 - x)H_ 2 ,o + 32/3(3 + 4/x - 3x - 4x 2 )(H li0 + H M ) 



+ ( 1 + x) (448/3H ,o + 64H 2 , + 64H 2 ,i )] } 
+ L 2 ,{c F (n / + l) 8/9(3 + 4/x -3x-4x 2 ) + 16/3(1 +x)H 



+C F 



1 6/3 (3x + 4x 2 )H - 4/3 (23 - 23x - 24^ 2 - 24^ 2 x) - 1 6/3 (3 + 4/x - 3x - 4x 2 )Hi 



■16(l+x)(Ho,o + 2H 2 ) 



+C A C F 



16(1 - 2x)H ,o - 8/3(11 + 4/x + 8x)H +4/3(45 - 44/x- 45x + 44x 2 + 12^ 2 + V2JC&) 



-8/3(3 + 4/x-3x-4x 2 )H! -16(1 +x)H 2 }} 



32/9(10 - 1/x - x - 8x 2 - 3^2 - %ix) - 16/9(3 + 15x - 4x 2 )H Q 



■16/9(3 + 4/x - 3x - 4x 2 )Hi - 32/3( 1 +x) (H 0;0 + H 2 ) 



+C F rif 



8/9(3 + 15x - 4x 2 )H - 16/9(10 - 1/x - x - 8x 2 - 3^ 2 - 3^ 2 x) 



-8/9(3 + 4/x -3x-4x 2 )Hi - 16/3(1 +x)(H , + H 2 ) 



+C F 



8/3(3 + 12x - 8x 2 )H 2 +4/3(33 - 33x - 18^ 3 - 18^ 3 x - 6^ 2 - 24£ 2 x + 16^ 2 x 2 ) 



57 



+2/3(41 - 53x - 64x 2 - 36^ 2 - 36^ 2 x)H + 2/3(103 - 8/x - 31x - 64x 2 )H! 
+ 16/3(3x-2x 2 )H , + 8/3(3 + 4/x-3x-4x 2 )(H 1 , + 2H u ) 

+(1 +x)(8H OA o + 24H 3 + 16H 2 , + 32H 2jl ) 

+C A C F 88/3(1 +4x)H 0j0 - 4/3(5 -40/x-41x + 76x 2 )Hi +8/3(14 + 4/x + 23x-4x 2 )H 2 

+4/3(15 - 23 /x + 195x - 187x 2 - 36^ 3 x - 8/x^ 2 - 28^ 2 - 46^ 2 x + 8^ 2 x 2 ) - 16(1 - 2x)H 0;0 ,o 
-4/3 (63 + 4/x - 1 05x + 44x 2 + 36^ 2 x)H + 8 /3 (3 + 4/x - 3x - 4x 2 ) (H i j0 + 2H u ) 

+48xH 3 + 16(1 +x)(H 2 , + 2H 2a ) 



L Q {c F (n f + 2) 



8/9(3 + 4/x-3x-4x 2 ) + 16/3(1 +x)H c 



+C F 



8 / 3 (3x + 4x 2 )H - 2/3 (23 - 23x - 24^ 2 - 24^ 2 x) - 8 / 3 (3 + 4/x - 3x - 4x 2 )H i 



-8(l+x)(H 0i0 + 2H 2 ; 



+c A c F 



16(1 - 2x)H ,o - 8/3(11 + 4/x + 8x)H +4/3(45 -44/x-45x + 44x 2 + 12^ 2 + 12^ 2 x) 



-8/3(3 + 4/x-3x-4x 2 )H! -16(1 +x)H 2 }} 



©10) 



Here c^' t "'(x,n f + 1) is the three-loop pure-singlet quark coefficient function at /u = Q as given 



,(3),ps, 
'2,q \ A "> n f 

in Eq. (B.10) of Ref. ifTTl . while «gJp S denotes the n^-independent part of the massive OME 

approximated in Eqs. (13.521 ) and (I3.531 ). Again, the terms proportional to flu in C2 3 j' ps (x,n^-+ 1) 
- cf. Eq. (4.12) in Ref. IfTTl - are omitted following the discussion in Section 1331 Expanding in 
powers of (X s (nf) instead of (X s (n f + 1) leads to the additional term 



r(3),PS, 
l 2,q 



m i ,v (x, Q\m 



2 2 a s {n f +\) -^a s (n f ) 



H 



g' ps (x,e 2 ,m 2 )-^L,//g' ps (x,e 2 ,m 2 ) 



dBii) 
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